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Predicting CT Contrast Enhancement With Feedback 

Fi.eld of the Invent:lon 

This invention relates generally to an apparatus 
and method for predicting organ specific contrast 
enhancement prior to computed tomography scanning of a 
patient. Specifically, this invention relates to a 
5 ccxnputer simulation of contrast agent transport 

throughout the body to predict organ specific enhancement 
in patients of variable height and weight subjected to 
various contrast injection protocols to enable operator 
selection of an appropriate injection protocol prior to 
10 commencing the scan and to use measurements of actual 
contrast agent transport throughout the body after 
injection as feedback to verify and calibrate the 
predictions. 

Backarotind of the Inv ention 

15 Computed tomography (CT) is a widespread 

diagnostic imaging method which measures the x-ray 
attenuation coefficient of matter. This x-ray 
attenuation coefficient is depicted in terms of 
5 Hounsefield Units (HU). During a CT scan, a collimated 
X-ray beam is directed on the patient and the attenuated 
remnant radiation is measured by a detector whose 
response is transmitted to a computer. The computer 
considers the location of the patient and the spatial 

10 relationship of the x-ray beam to the region of Interest. 
The computer analyzes the signal from the detector so 
that a visual image can be reconstructed and displayed on 
a monitor. The image can then be viewed or stored for 
later evaluation. 

25 Hounsefield Units reflect the relative absorption 

of CT x-rays by matter, the absorption being related to 
the atomic number, electron density, physical thickness 
of that matter, and the energy spectrum of the x-rays. 
Because of the similarity in electron density of various 
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tissues in the body, CT scans sometimes result in poor 
imaging. In an attempt to obtain better results in such 
circumstances, a contrast agent, such as iodine, can be 
injected in the patent's blood stream to change the 
5 relative radio-density of the tissues, and improve 
overall diagnostic efficacy. 

When using a contrast agent, it is extremely 
important to coordinate the time of the scan with the 
time of greatest levels of contrast in the region of 
10 interest, in some instances with respect to a threshold 
value. Because the contrast agent is injected into the 
blood stream, many physiological factors can affect the 
start time and duration of a sufficient level of contrast 
in the region of interest. For example, because the 
15 cardiovascular system provides the means for circulation 
of contrast agent throughout the body after it is 
injected into the blood stream, a patient's cardiac 
output can have a significant effect on the distribution 
of the contrast agent as well as the time taken for the 
20 contrast agent to reach a particular organ or vessel. 

Current understanding of intravenous contrast 
enhancement is further complicated by multiple 
interacting factors including contrast agent type, volume 
and concentration, injection technique, catheter size and 
25 site, scanning technique, patient characteristics and 
tissue characteristics. Of these factors, all of which 
have influence on contrast enhancement, the variables 
which cannot be controlled are those related to the 
patient. These include age, gender, weight, height, 
30 cardiovascular status, renal function and other disease 
status. m the past ten years, many clinical studies 
testing various intravascular contrast agent and 
injection protocols have been reported. However, in many 
respects, contrast enhancement still relies heavily on 
35 the experience and intuition of the physician rather than 
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rigorous, quantitative analysis of the mechanism of 
contrast enhancement. 
gi.inni»*T- y o£ tb ? Invention 

The present invention provides a method of and an 
apparatus for predicting tissue specific CT contrast 
enhancement in a patient with a specific body habitus 
subjected to different contrast injection protocols. The 
5 method is preferably implemented in a computer program, 
and the computer itself may be used to also control the 
CT scan in accordance with an operator's choice. Such a 
physiological model of contrast enhancement has many 
potential clinical applications. 
j^Q The present Invention utilizes a compartmental 

model of the human cardiovascular system and assigns 
differential equations describing mass transport to each 
compartment of the model. Regional circulation 
parameters such as blood volume, regional blood flow and 
15 extracellular fluid volume were estimated using available 
data to provide input to the equations. Local tissue 
structures such as organs and vessels were modeled 
mathematically to describe the distribution and 
dispersion of an intravascularly-administrated contrast 
20 agent. A global model was then formed by integrating the 
regional circulation parameters with the models of local 
tissue structures. 

The present invention, which is preferably 
implemented in a computer program, allows an accurate 
25 prediction of the time varying distribution and 

concentration of contrast in the body. This in turn 
allows an operator to predict the time and duration of 
maximum enhancement in a specific organ or tissue in a 
patient for a particular injection protocol. Most 
30 importantly, the operator can use the present invention 
to predict the time a scan should be started and the 
duration of the scan based on output data from the 
program. This output can take the form of a data stream 
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or can be a graph of contrast enhancement, versus time. 
With more advanced generations of CT machines, such as 
spiral and helical CT scanners as are known in the art, a 
typical scanning procedure can be completed within 
5 approximately 30 seconds* The present Invention enables 
an operator to choose an Injection protocol to ensure 
that the entire scan takes place during a period of 
maximum enhancement, and while enhancement exceeds a 
suitable threshold. 
10 In the prior art, there are devices which monitor 

and output contrast enhancement levels for a region of 
Interest. Using these prior art devices, an operator 
Injects the contrast agent Into the patient, views the 
output, and determines when to begin a scan based on when 
15 the enhancement level attained In the patient ' s region of 
Interest becomes acceptable* The prior art devices 
require the Injection of the contrast agent and low dose 
x-rays of the region of Interest. For example. In the 
prior art, the Injection of a contrast agent Is started 
20 and the prior art device monitors at regular Intervals 
the enhancement level In the region of Interest and 
provides an output. The operator can view the output and 
then decide when to begin the scan. 

Prior art techniques to determine the enhancement 
25 level in an organ Include low dose pre- scanning of the 

organ during injection of the contrast. One such example 
is U.S. Patent No. 5,459,769 to Brown. In Brown, after a 
short delay upon initiating the injection, low-dose x- 
rays are taken of the organ to be scanned, images are 
30 constructed from the low-dose x-rays and the Images are 
displayed for the operator to determine when the 
enhancement level in the organ is sufficient to begin a 
full dose scan. 

The present Invention is a significant improvement 
35 over the prior art in that it allows prediction of 
contrast enhancement levels and duration of those 
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enhancement levels prior to injection of the contrast 
agent and without the need of low dose x-rays. Moreover, 
because different injection variables, such as rate and 
concentration, can alter the enhancement levels, the 
5 present invention allows calculation of various 

alternatives to choose the best injection scheme for a 
particular patient. The prior art devices do not assist 
the operator in determining whether a particular patient 
having a specific body habitus will even acquire a 

10 desired threshold level of enhancement with a given 
injection protocol. 

Therefore, if the threshold level of enhancement 
is never reached because of the patient's specific body 
habitus or an improper injection protocol for the 

15 particular patient, a scan cannot be completed and the 
entire process must be repeated, including a second 
injection. Even then, the operator cannot be certain 
that the revised injection will ever reach a desired 
threshold level of enhancement based on that patient's 

20 specific body habitus and injection protocol, nor whether 
the threshold level, if attained, will be maintained 
during the entirety of a desired scan duration. This is 
particularly important for scans of certain tissues whose 
contrast enhancement behavior is complex, as explained in 

25 greater detail below. 

The present invention also allows an operator to 
adjust the collimation or slice thickness and CT table 
speed to optimize a scan. During a CT scan, a patient 
lies on a table which moves through the CT scanner from 

30 head to toe vertically, and over the selected region of 
interest. The collimation or slice thickness is the 
thickness of the slice of the patient's body that is 
transaxially scanned. The table can usually be moved at 
a rate per second of up to two times the collimation 

35 thickness. Using the method and apparatus of the present 
invention, an operator can optimize the collimation rate 
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and table speed. For example. If -there is a limited 
period of threshold enhancement, the operator knows that 
an increased table speed or an Increased collimation 
thickness must be used to ensure that the entire scan is 
5 completed during the time period of maximum enhancement. 
Customizing the scan is less precise with the prior art. 

In the present invention, the computer program 
allows an operator to determine how long the predicted 
enhancement level exceeds the threshold enhancement. 

10 That information provides the operator with a means to 
adjust various scan parameters such as scan duration, 
scan start time and table speed, to be certain that the 
scan takes place during a period when the predicted 
enhancement is above a threshold enhancement. Thus, the 

15 invention was a significant improvement over the prior 
art* 

Building on the invention as described above, the 
inventors herein have improved it by providing for a 
method of using the predicted enhancement levels to 
20 optimize injection protocol and adjust the injection 
parameters to increase or decrease the time that the 
predicted enhancement level exceeds the threshold. The 
inventors have also improved it by providing a means for 
determining the optimum scan start time when the 
25 predicted enhancement level exceeds the threshold for a 
period much greater than the scan duration. 

The present invention can be implemented in many 
ways including a separate computer or integrated with the 
computer of a CT machine. All that is required is a 
30 computer having the invention programmed therein 

integrated with the controls of the machine. The present 
invention can also be implemented in a contrast injector 
system which is equipped with a computer for predicting 
contrast enhancement for given inputs of patient 
35 parameters and contrast injection protocol. In this way. 
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adjustmen-cs to the Injection protocol are readily made 
using the Injector. 

The contrast Injector having a computer with the 
invention programmed therein can be operated as a fully 
5 integrated system with a CT scanner or as an Independent 
system* When the present Invention is used as a part of 
a computer system integrated with both an injector system 
and a CT scanner equipped with low-dose pre- scan CT, the 
optimal set of scan parameters can be adjusted based on 
10 actual enhancement measurements acquired with low-dose 
scanning . 

The present Invention is capable of using standard 
values for variables which influence enhancement levels 
and also allow input of patient specific values. For 

15 example, a particular patient habitus may be such that 
the standard values for variables such as blood volume, 
blood flow etc., will not provide an accurate prediction 
of enhancement levels. The invention utilizes several 
methods to resolve such situations. One method provides 

20 for the input of patient specific information to 

customize the operation to the particular patient. This 
includes patient specific variables such as weight, age, 
height and gender. These variables can be measured and 
input to adjust the standard variables accordingly. 

25 On occasion, other variables which are not readily 

measurable may need to be modified. As is well known in 
the art, cardiac output cannot be measured as readily as 
height or weight. Of course, a patient with a known 
history of heart failure or increased age will most 

30 certainly have a cardiac output below normal. If this is 
the case, the invention of the parent allows adjustment 
of the standard variables accordingly. 

Another aspect of the method allows the operator 
to choose several alternative values for cardiac output 

35 and generate a family of predicted enhancement curves for 
each value. After injection of contrast agent is 
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s'tar'ted, actual measurements of enhancement can be 
compared to the Initial portions of the family of curves 
to determine which family member most closely resembles 
the actual results. In this way, early in the scan and 
5 before the threshold has been reached, a choice of which 
curve to utilize to best predict when the scan should 
occur can be made. This choice can be made by the 
operator or automatically by the computer. 

The inventors herein have improved upon this 
10 method by using predicted aortic enhancement levels 
compared with sequential measurement of actual aortic 
enhancement levels using low-dose pre-scanning as an 
indicator of unknown patient specific parameters such as 
the cardiac output of the patient. 
15 In the present invention, the predicted 

enhancement levels are computed for a given set of 
patient specific parameters in an injection protocol. 
Prior to any injection of the contrast into the patient, 
the operator enters the patient specific parameters and 
20 the injection protocol into the computer which can be 
integral with an injector, a CT scanner or be a stand- 
alone personal computer. The method of the present 
invention provides an output from the computer which 
gives the predicted tissue enhancement level of the 
25 tissue to be scanned as a function of time. Based on 
that output, the operator uses the present invention to 
modify the injection protocol in order to ensure that the 
predicted tissue enhancement function exceeds a threshold 
level for at least as long as the desired scan duration. 
3Q If the predicted tissue enhancement function is 

shown by the output to be substantially greater than the 
threshold level, the present invention allows the 
operator to modify the injection protocol to decrease the 
volume or flow rate. If the predicted enhancement 
35 function does not meet or exceed the threshold level or 
does not exceed the threshold level for the length of the 
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scan duration, the present invention allows modification 
of the injection protocol by increasing the flow rate or 
volume of contrast. The program iterates to again 
compute and provide a revised predicted tissue 

5 enhancement function output based on the revised 
injection and/or rate. 

If the predicted tissue enhancement function is 
known to be above the threshold level for a period of 
time substantially greater than the desired scan 
10 duration, the operator can use the method above to 

decrease the Injection rate or volume of contrast or both 
to reduce the predicted enhancement level to a range 
which still satisfies the scan parameters, thereby saving 
contrast agent and minimizing any potential side effects 

15 to the patient. in the alternative the method provides 
for prediction of the optimal temporal window for 
performing the scan within the period that the 
enhancement exceeds the threshold. 

After the operator is satisfied that the injection 

20 protocol chosen and the patient specific parameters will 
produce an acceptable enhancement level as shown by the 
output of the present invention, the operator can further 
increase the accuracy of the prediction by predicting and 
low-dose monitoring enhancement in a region of Interest 

25 and comparing the prediction with the actual measurement 
from the monitoring to update or revise the predicted 
tissue enhancement function. in this way, the invention 
uses feedback from actual enhancement measurements to 
fine tune the predictions. 

30 To practice this aspect of the invention, the 

operator performs a base line scan over a distinct region 
of interest. The base line scan is a low-dose or partial 
scan in which the x-ray dosage is reduced substantially 
less than a typical scan. After completing the base line 

35 scan, the operator begins the injection of the contrast 
agent. After initiating the injection, the operator 
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performs low-dose pre-scan of the region of interest, 
such as the aorta, to obtain actual aortic enhancement 
levels. The pre-scan is virtually Identical to the base 
line scan described above and is also a low-dose scan and 
5 can use less than a full revolution of the gantry of a CT 
scanner. Using the actual measurement of enhancement in 
the aorta, the present invention can calibrate the model 
for patient parameters such as cardiac output and provide 
revised tissue enhancement predictions. 

10 The present invention has particular application 

where the organ or vessel being scanned is incapable of 
maintaining threshold enhancement levels for a sustained 
period. One such example is CT angiography. In CT 
angiography, a CT scan is taken of a blood vessel or 

15 vessels. Unlike organs, blood vessels do not maintain 
high enhancement levels over time and the timing of the 
scan is critical. CT angiography is performed in the 
prior art by injecting a test dose of contrast agent and 
measuring with low dose x-rays the elapsed time for the 

20 contrast agent to reach the region of interest. 

Thereafter, a full dose of contrast agent is injected and 
a scan is initiated after lapse of the previously 
measured time delay. However, there is no guarantee of a 
particular enhancement level being attained or sustained 

25 as required to achieve a successful scan. Using the 

invention comprising the method and apparatus disclosed 
herein, one can more accurately predict not only the time 
delay, but also the degree of enhancement and its 
duration • 

30 While the principal advantages and features of the 

invention have been described above, a greater 
understanding of the invention may be attained by 
referring to the drawings and the description of the 
preferred embodiment which follow. 
BT-i ^f Description of the Drawings 
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Figure 1 is a diagram showing the components of a 
complete CT scanner system and a computer control 
console; 

Figure 2 is a schematic diagram of the major 
5 organs of a human cardiovascular circulation systems- 
Figure 3(a) is a block diagram of a single well 
stirred compartment with an input having a constant input 
concentration and input flow rate Qj, the compartment 
having a volume V and an output with a concentration C, 
10 and an output flow rate Q^; 

Figure 3(b) is a graph of the input concentration 
of the input in Figure 3(a) and a graph of the 
corresponding output concentration C„ over time; 

Figure 4(a) is a bloclc diagram of an organ modeled 
15 in three spaces: intravascular (IV), extracellular (EC), 
and intracellular (IC); 

Figure 4(b) is a block diagram of the IV and EC 
spaces of Figure 4(a) detailing the mass transfer rate 
( dM/dt ) therebetween ; 
20 Figure 5 is a block diagram of the global 

cardiovascular model of the body; 

Figures 6 is a flow chart showing the method steps 
for determining predicted contrast enhancement level; 

Figure 7 is a flow chart of a subroutine of the 
25 method of Figure 6 for operator designation of patient 
information and contrast protocol information; 

Figure 8 is a flow chart of a subroutine of the 
method of Figure 6 which assigns a differential equation 
to each element of the cardiovascular model; 

Figure 9 is a graph showing the linear 
relationship between enhancement in Hounsefield Units 
(H.U. ) and concentration of Iodine (I mg/ml); 

Figure 10 is a graph showing simulated (10b) and 
empiric (10a) aortic and hepatic enhancement using the 
35 biphasic-low flow rate injection protocol given in Table 
6; 



30 



wo 97/12550 



PCTAJS96/15680 



12 

Figure 11 Is a pair of graphs showing simulated 
(lib) and empiric (11a) aortic and hepatic enhancement 
using the unlphaslc-low flow rate Injection protocol 
given In Table 6; 
5 Figure 12 Is a graph showing simulated (12b) and 

empiric (12a) aortic and hepatic enhancement using the 
unlphaslc-hlgh flow rate injection protocol given in 
Table 6; 

Figure 13 is a graph showing simulated aortic and 
10 hepatic enhancement curves generated by the Invention for 
hypothetical patients weighing 110, 150, 200 and 250 
pounds; 

Figure 14 Is a graph showing simulated aortic 
(14a) and hepatic (14b) enhancement curves for 
15 permeability (PS) values of -1, 1.0, 20 and infinity; 

Figure 15 is a data stream output generated by the 
invention showing predicted aortic and hepatic 
enhancement levels in a hypothetical patient with 
standard blood volume and standard cardiac output using 
20 unlphaslc-hlgh flow rate injection protocol In Table Six; 
Figure 16 Is a graph output generated by the 
invention showing predicted aortic and hepatic 
enhancement levels versus time using the data in Figure 
15; 

25 Figure 17 Is a drawing of a power injector system 

with Its various components; 

Figure 18 is a graph output generated by the 
present Invention showing the predicted enhancement 
function versus time superimposed over a line 

30 representing a threshold enhancement level; 

Figure 19 Is a graph output generated by the 
present invention showing the predicted enhancement 
function and Identifying those intervals for which the 
predicted enhancement function exceeds the threshold 

35 value; 
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Figure 20 is a table showing an output of the 
present Invention using the process steps In the flow 
chart of Figure 21; 

Figure 21 Is a flow chart showing the method steps 
5 of the present Invention used to calculate an Injection 
flow rate and volume which provide a predicted 
enhancement function which exceeds a threshold value for 
a period greater than the scan duration and the method 
steps to select an optimum scan interval within the 
10 period; 

Figures 22A and 22B are graphs generated using the 
present invention showing four predicted aortic 
enhancement curves 22A and four different predicted 
hepatic enhancement curves 22B corresponding to four 

15 alternative cardiac outputs; 

Figure 23 is a graph showing predicted and actual 
aortic enhancement along with predicted and actual 
hepatic enhancement for a given patient; 

Figure 24 is a table generated using the present 

20 Invention showing four predicted aortic enhancement 
levels for four alternative cardiac outputs and the 
actual predicted enhancement levels at specified elapsed 
times after injection with the corresponding area under 
the enhancement curve (AUC) calculations; 

25 Figure 25 is a table generated using the present 

invention showing predicted hepatic enhancement levels at 
specified elapsed times after injection using the cardiac 
output calculated from the table of Figure 24; and 

Figure 26 is a flow chart showing the method steps 

30 of the present invention to select the optimum onset of 
hepatic scanning from a table of predicted optimum onset 
times correlated to actual measured times for the actual 
enhancement curve to achieve an AUC of a predetermined 
amount . 

n^<-, *<lftd Description o f the Preferred Embodiment 
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Computed Tomographic (CT) scanning is an 
invaluable radiologic diagnostic tool. The major 
components of a conventional CT scanner are shown in 
Figure 1 . The CT scanner 10 contains the x*ray tube and 
5 detector array. Power is supplied by a high voltage 
generator 12 controlled by scanner electronics 14 and 
scanner service module 16- The patient support and 
positioning couch 18 is moveable to transport the patient 
through the scanner 10. The scanner 10 and voltage 
10 generator 12 receive electronic commands from the 

operating console 20 and transmit data to the computer 
system 22 for image production and analysis. The 
operating console 20 usually contains an interactive 
keyboard 24 and CRT monitor 26. 
15 Many radiographic procedures, including CT scans, 

require an injection of contrast medium under specific 
control conditions. For example, CT scanning requires a 
high degree of control over the injection of the contrast 
agent and the parameters of the injection protocol in 
20 order to maximize the accuracy of the scan. It would be 
difficult to consistently perform these injections by 
hand. Therefore, these injections are usually performed 
with a mechanical -type device known in the art as power 
injectors or injector systems- Injector systems allow 
25 significant control of the injection of contrast into a 
patient prior to a CT scan. 

Injection systems for CT scans have several basic 
components. An example of an Injector system Is shown 
generally as 128 in Figure 17. The injector system 
30 includes a control unit 130, mounted on a pedestal 132, 
an injector head 134 and arm 136 integral to the unit. 
The entire injector system 128 can be wheeled about as a 
unit for use with the CT scan machine as shown in Figure 
1. The control unit of the Injector comprises a control 
35 panel for setting up the injection and a display for 
displaying instructions and data. The controls and 
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Indicators present on the control panel will vary with 
the type of options available on the system. One example 
of an injector system is the Mark V Plus injection system 
manufactured by Medrad of Pittsburgh, Pennsylvania. 
5 The injector head accommodates a syringe and 

provides for power injecting. The arm connects the head 
to the console permitting easy movement for loading or 
injecting. The height allows placement of the head over 
the patient during a CT scan. In the alternative, the 

10 injector head and arm can be mounted to the CT scan 

table, or overhead from a ceiling mount, or wall mounted 
near a CT system with the control panel being located on 
a console integrated with the controls of the CT scan. 

The injector head can be either floor mounted or 

15 track mounted on an overhead system. In most cases, the 
control panel is mounted inside the CT scan control booth 
for operator safety. Most injector systems also contain 
a microprocessor and memory for storing computer programs 
which can be recalled when needed. A warming system may 

20 also be included in an injector system which warms and 
maintains the contrast medium at or near body 
temperature. This will help reduce viscosity of contrast 
medium, resulting in a decrease in resistance of contrast 
medium flow and decrease in a patient discomfort 

25 experienced during injection. Most injector systems use 
a power high pressure mechanism, such as an electric 
drive motor coupled to a jack screw, to drive the piston 
in or out of the syringe and deliver the contrast. 

Injector systems can also be interfaced with a CT 

30 imaging system. This interfacing of an injector to a CT 
imaging system allows variations such as causing the 
injector to be triggered by the imaging system or causing 
the imaging system to be triggered by the injector. Bi- 
directional controls are also available to allow either 

35 device to control the other and allow an operator to 

choose how to sequence and time the devices when they are 
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connected. For example, the CT scan Imaging system can 
have a control which allows It to trigger the Injector 
and the Injector can have a control which allows the 
Injector to send a signal to the CT Imaging system to 
5 trigger the start of a scan* The details of Interfacing 
the injector system with an imaging system and the 
controls available are well-known in the art. 

The control panel on an injector system is used to 
set the parameters of the injection sequence. It usually 
10 consists of an alpha-numeric keyboard and buttons for 
various input parameters as well as several display 
windows Including a window for displaying system 
messages. The control panel accepts and displays 
injection parameters, displays Injection results and 
15 other messages related to the control of the injector 

system. The control panel allows the operator to program 
the injector system to control various parameters of the 
injection process including the flow rate, volume, 
injection duration, injection pressure, and injection 
20 delay. The flow rate is defined as the delivery rate of 
the contrast (amount delivered per unit of time). The 
flow rate is dependent on the viscosity of the contrast 
agent, the length and diameter of the catheter, and the 
Injection pressure. The particular flow rate chosen for 
25 a specific procedure is governed by the procedure itself, 
the vessel entered and the patient habitus. Flow rates 
can vary from as low as •! ml/s to as high as 40 ml/s, 
depending on these factors. 

The present invention Is preferably implemented in 
30 a computer program. Because most CT scanners and 
Injector systems utilize computers to control their 
operation, the present invention could be easily 
Integrated therewith. In that fashion, the operator 
could run the program prior to the injection or prior to 
35 the scan and the computer can determine the optimum 

Injection and scan parameters and using the determined 
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values complete the injection and the scan accordingly. 
In this way, the injection system computer or the CT scan 
computer (or the shared cwnputer on a combined syst»m) 
can determine as well as implement the injection protocol 
5 and the scan parameters. In the alternative, a separate 
computer which contains the program could be utilized. 
The scan parameters and injection protocol could be 
determined by running the program in the separate 
computer and then input into the injector system computer 
10 or CT scan computer by the operator or through computer 
data transfer methods. 

The invention utilizes a model of the human 
cardiovascular system to describe mass transport of 
contrast agent throughout the body. The cardiovascular 
15 system provides the means for circulation of contrast 
agent throughout the body after it is injected into the 
bloodstream. The human cardiovascular system is very 
complex and has numerous controlling mechanisms, 
including neuronal, hormonal and psychological controls. 
20 A simplified human cardiovascular system as shown 

schematically in Figure 2 consists of the heart, vascular 
networks, and key organs which serve as reservoirs. 
Normal blood volume and flow distribution throughout the 
body are well established in the prior art and are given 
25 in the Tables 1 and 2 (All Tables are shown in Exhibit A 
attached hereto and incorporated herein by reference). 

Based on well known information, the model assumed 
that the average blood volume was 5 liters. This 
includes 3 liters of plasma and 2 liters of red blood 
30 cells. The average cardiac output was also estimated 
from known sources to be 6.5 liters per minute. These 
values were used to describe a standard model of the 
cardiovascular system. However, the method and apparatus 
of the invention allow these values to be adjusted 
35 according to the patient's age, gender, weight and height 
using standard nomograms outlined below. 
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Because conlrras^t agent diffuses passively from the 
bloodstream across the capillary membrane Into the 
extravascular space, the distribution of fluid throughout 
the body was Included in the cardiovascular model* The 
5 amount of total body water (TBW) in an adult of average 
weight (70 Kg) was assumed to be 40L. TBW was divided 
into two major components. Intracellular fluid (ICF) and 
extracellular fluid (ECF). The ECF was further divided 
into several smaller compartments including interstitial 

10 fluid, plasma, and cerebrospinal fluid. The interstitial 
fluid is the largest compartment and lies in the 
lymphatics and the spaces between cells. 

The ECF volume is usually estimated with dilution 
methods in which a substance is injected into the blood 

15 stream and diffuses throughout the entire extracellular 
fluid compartment with little entering into the cells. 
However, an ideal substance for such dilution studies has 
not been identified, and measurements for a 70kg adult 
have ranged from 9L to 22L depending on the substances 

20 used. The size of the measured ECF decreases with 

increases in the molecular weight of the substance used. 
The apparent volume distribution of iohexol has been 
reported to be .27 1/kg. Thus, for a 70Kg adult, this 
equates to an ECF of 18. 9L. In the model this value for 

25 ECF volume was used which includes a plasma volume of 
3.0L. The overall estimated distribution of body fluid 
used in the cardiovascular model is summarized in Table 
3. 

The detailed distribution of fluid In a local 
30 organ was estimated from the standard mass of an organ 
and its water content. The volume of the total systemic 
capillary bed is estimated to be about 300 ml. However, 
a detailed breakdown of capillary volumes in different 
regions is not available. In addition, the number of 
35 capillaries within an organ varies considerably from one 
organ to another. It is believed that the regional 
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10 



capillary volume is directly proportional to a regional 
blood flow and the cardiovascular model applied this 
assumption. These values are likely overestimated in 
highly perfused organs such as the kidney and the liver 
but this did not hinder the performance of the model. 
Table 4 shows the regional capillary volumes in the 
systemic circulation estimated from the regional blood 
flow values given in Table 2. 

A calculation of the regional distribution of the 
extracellular and intracellular fluid was also necessary 
for the invention. The regional distribution of total 
body fluid can be calculated from the known mass of each 
organ and its water content, assuming a density of 1.0 
g/Ml- The weight and percent of water content of the 
15 visceral organs are shown in Table 5 along with their 
total fluid value minus the capillary volume. Without 
available information, the model assumes 70% water 
content in the stomach, spleen and intestine. The lung 
consists of 50% parenchyma and 50% non-parenchyma tissues 
20 whose capillary volumes are 150 ml and 5 ml, 
respectively . 

The total body fluid of the upper extremities, 
trunk, and lower extremities was calculated by 
subtracting from the total body water volume (40L), the 
25 blood volume (5L), and the total fluid of the visceral 
organs minus the capillaries (4,726 mL). The mass ratio 
of the lower extremities and trunk to the upper 
extremities is about 4:1. Thus, the total body fluid of 
the upper extremities becomes 6,055 mL and that of the 
30 lower extremities and trunk, 24,219 mL. 

Table 3 shows the overall volumes of ICF fluid and 
ECF fluid as 19. IL and 15.9 L. respectively. However, 
regional distribution of the ICF and ECF is not shown. 
Some tissues such as the skin, adipose tissue, G-l tract, 
35 and liver have larger extracellular to intracellular 
fluid ratios than, for example, muscle. As no data 
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regarding such fluid ratios are available. It was assumed 
the ratio o£ BCF to ICF to be the same In all body 
regions. For example, the ECF and ICF volumes of the 
liver were estimated as 524 and 629 mL respectively. 
5 After regional blood flow, blood volume, and 

distribution of body fluid were estimated, local 
structures were modeled mathematically to describe the 
distribution and dispersion of intravascularly 
administered lodinated contrast agent within local 
10 regions. The blood vessels are viscoelastic with complex 
mechanical properties to accommodate pulsatile blood flow 
and various pressure gradients. Although the blood flow 
in large vessels is generally streamlined, some mixing 
occurs within the blood vessel because of molecular 
15 diffusion, flow pulsabillty and convections at multiple 
branching points. The dispersion may be even greater in 
smaller and low pressure vessels. To simplify the model, 
blood vessels were represented as rigid structures 
without directly incorporating their dynamic pulsatile 
20 properties in the. 

A blood vessel could be analyzed in the 
cardiovascular model as a simple conduit without any 
longitudinal mixing. This is known in the art as "plug 
flow." In this type of model, each artery and vein is 
25 divided into segments, and blood enters as plugs for each 
heartbeat and displaces an equal volume of blood without 
any longitudinal mixing. The major problem with this 
approach is excessive demands on computer memory required 
to store the history of each segment throughout the 
30 circulation. An alternative approach is to consider a 
blood vessel as a well-stirred compartment or well-mixed 
pool of blood. This approach simplifies computation, 
requires far less computer storage and has been shown to 
perform as well as the plug flow model in the prior art. 
35 Thus, in the cardiovascular model, the heart and blood 
vessels were analyzed as well-stirred compartments. 
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A single, well-stirred compartment contains a 
constant volume, V, with a single inlet flow and a single 
outlet flow as shown in Fig. 3(a). 0^ and represent 
the input and output volumetric flow rates of the blood, 
5 respectively. The input and output flow rates are the 
same in a constant volume compartment (Q-Qi-Qo). 
represent the input and output concentrations of contrast 
agent, respectively. Since we assume the compartment to 
be well mixed, the concentration within the compartment 
10 is the same as that of the output. A mass balance of the 
concentration is described by Pick's Principle, shown 
schematically in Figures 3(a) and 3(b) in the following 
equation : 

V*dCo/dt-Q( Ci-C, ) . 
15 For a given volume, V, a given volumetric flow 

rate Q, and a given input concentration, C^, we can 
estimate the output concentration, C„, by solving this 
differential equation. The net effect of a well mixed 
compartment is to disperse the input concentration over 

20 the compartment resulting in more broadly distributed 

output concentration over time. For constant flow rate, 
Q over a fixed time interval, T, the input concentration 
given as a step function is mathematically transformed to 
the output concentration curve as shown in Fig. 3(b). 

25 The transformation is described mathematically as two 
exponential functions of V, 0 and T. The output 
concentration curve is broader temporally than the input 
concentration curve, and a central peak is present. 

Modeling an organ is more complex than modeling a 

30 blood vessel because the contrast agent is no longer 

confined in the intravascular space and permeates through 
the capillary membrane into the extravascular space. The 
simplest approach to modeling an organ is to assume that 
it also is a well-stirred compartment. However, the 

35 single compartment organ model does not address 

differences in the exchange of contrast agent along 
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subcompartments within an organ and is llml-ted in 
describing the behavior of substances with different 
transcapillary permeabilities. A common alternative 
approach used in the prior art to investigate the 
5 distribution of chemotherapeutic agents throughout the 
body involves splitting each organ into three well known 
spaces: the capillary or intravascular space (IV)^ the 
extracellular space (EC), and the intracellular space 
(IC). This is shown schematically in Fig. 4(a). For a 

10 given organ, each of these three spaces was modeled as a 
single, well-mixed compartment. Diffusion through 
membranes, either active or passive, permits exchange of 
substances along the spaces within the organ. However, 
because iodinated contrast agent does not penetrate into 

15 the cells, only the intravascular (IV) and extracellular 
(EC) compartments were considered and the intracellular 
(IC) compartment was ignored. 

Transcapillary exchange of substances between the 
intravascular and extracellular compartments can be 

20 described by Fick's Law of Diffusion and is shown 
schematically in Figures 4(a) and 4(b). The mass 
transfer rate (dM/dt) is proportional to the diffusion 
coefficient (D), the surface area (S), and the 
concentration difference (C^-C^) for a given membrane 

25 thickness (dX) as represented by the following equation: 

dM/dt=DS( C^-C^ ) /dX - 
For a thin membrane, the mass transfer rate is simpler 
such that permeability (P) is commonly used to combine D 
and dX as a unit resulting in the following equation: 

30 dM/dt»PS(Ci-CJ 

To complete the mathematical model, two governing 
differential equations were applied to each organ. One 
for the intravascular space and the other for the 
extracellular space. The intravascular space had two 

35 transport components. The first component was obtained 
from blood flow related mass balance, i.e., the inflow of 
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contrast: agent minus the outflow. The second component 
was obtained from the mass balance related to the 
transcapillary exchange within the extracellular space. 
For the extracellular space, only one transport component 
5 was considered: the mass balance related to 

transcapillary exchange with the Intravascular space. 
These equations are as follows: 

Vi^*dCi,/dt-q( C, -Ci, ) -PS ( C^.-C.^ ) 
V.e*dC.,/dt=PS( Ci,-C., ) 
j_0 The global model, shown schematically in Figure 5, 

was formed by integrating the regional circulation 
parameters with the models of local regions. In the 
model, contrast agent was assumed to be Injected through 
an antecubltal vein, mixed In the right heart, 
15 distributed throughout the body and excreted by the 
kidneys according to the glomerular filtration rate. 

The residence time of contrast agent In an organ 
was estimated by the time duration of the contrast agent 
in the capillaries and ECF spaces. The residence time 
20 depends on the size of these spaces as well as the 
transcapillary exchange rate. When a substance Is 
confined to a blood vessel, the circulation time is 
measured by injecting rapidly a dye or radioactive tracer 
Into a peripheral vein and detecting the moment when it 
25 arrives at a sampling site. The volume of a blood vessel 
travelled by a substance is calculated by multiplying the 
volumetric flow rate and the circulation time. The mean 
circulation time from the antecubltal vein to the right 
atrium is approximately 6.9 seconds in an average adult. 
30 The time can range from 3 to 14 seconds. This is the 

temporal difference between the antecubltal and the right 

atrial injections* 

Intravascular contrast agents are eliminated from 
the body mainly by the kidneys. The process is rapid 
35 with approximately 50% of injected contrast agent being 
excreted within two hours presuming normal renal 
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function. The totial excretion rat:e of contrast agent Is 
obtained by multiplying the plasma concentration with a 
glomerular filtration rate, usually about 19% of renal 
plasma flow. Peak renal excretion Is closely related to 
5 peak plasma concentration, because renal plasma flow is 
relatively constant. 

Regional blood flow Is expressed according to the 
magnitude and direction of the flow. For example, the 
cardiac output is 6500 mL/mln., directed away from the 
10 right heart. In Figure 5, the right and left heart are 
represented by boxes by denoting well stirred 
compartments. Each blood vessel is represented by a 
circle surrounding a number which represents its volume 
in milliliters. Large blood vessels are further divided 
15 into multiple smaller compartments in series, typically 
20 mL for arteries, and 100 mL for veins: the voLxmae of 
systemic veins is about 4 to 5 times that of associated 
arteries. This division scheme in large vessels is 
rather arbitrary and was based on computational 
20 convenience. However, the total blood volume in a given 
blood vessel closely followed known physiological values. 

In Figure 5, each organ is sho%m as a box split 
into two sub-compartments, the upper number denoting 
Intravascular (capillary) volume and the lower number 
25 denoting extracellular fluid volume. The concentration 
of contrast agent in an organ is determined by the ratio 
of the total mass to the total volume of contrast agent 
within that organ. The total mass of contrast agent 
within an organ is calculated by summing the products of 
30 the concentrations and the volume in the Intravascular 
and extracellular spaces. The organ volume is obtained 
by adding the Intravascular (IV), extracellular (EC) and 
Intracellular (IC) volumes. 

A total of 104 ordinary differential equations 
35 were used to describe the cardiovascular model. These 
equations were solved using the numerical integration 
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programs of the fifth order Runge-Kutta method on a 
personal computer. Using a power Macintosh or IBM PC the 
computation took a few seconds to compute. The contrast 
concentration curve over time was calculated for each 

5 region by solving differential equations of the model for 
a given contrast Injection protocol and a hypothetical 
patient with variable weight, height, gender and age. 

Referring to Figure 6, the method Is shown as a 
flow chart. The first step In the method Is to call the 

0 patient /contrast information subroutine shown in Figure 
7. This subroutine accepts operator input of patient and 
contrast information which will affect physiological 
parameters of contrast enhancement • First, a 
permeability factor with a range from 2 to 10 is input. 

5 Guidance for selection of an appropriate permeability 
factor is given, infra. However, the Inventor has found 
that acceptable results are achieved upon operator 
selection of any number between 2 and 10. Next, a body 
size option Is input. At this point, the user has a 

to choice to use a standard model which will include a 5,000 
xnilllliter blood volume and standard cardiac output of 
6500 ml/min or a user may input specific information. If 
specific information is input, the standard blood volume 
and standard cardiac output are adjusted to conform to 

15 the patient's specific Information. Blood volume (BV) 
and cardiac output (CO) can be predicted from the weight 
(W) in pounds and height (H) in inches of a patient using 
regression formulae available in standard cardiovascular 
physiology references. The formula for an adult male 

)0 with a weight (W) ranging from 100 to 310 pounds, and a 
height (H) ranging from 60 to 74 inches is: 

BV - 33.164 * H°^" * W°*" - 1229. 
For an adult female with weight (W) ranging from 80 to 
290 pounds and height (H) ranging from 60 to 74 Inches 

}5 the formula is: 

BV - 34.85 * H®"^ * W^*25 . ^954^ 
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For an adult male or female, the cardiac output 
(CO) Is given by the formula: 

CO - 36.36 * H"'" * W**" 
In the model, an adjustment to these variables was 
5 made as follows. The ratio of the predicted blood volume 
to the standard blood volume was calculated. This ratio 
was then applied to the regional blood volxime and 
extravascular fluid volume in the cardiovascular model so 
the entire body fluid volume was corrected. The cardiac 
10 output and regional blood flow were also modified in the 
same fashion. Consequently, the regional blood flow, 
blood volume, and distribution of body fluid in the model 
can be adjusted for subjects of different body weight, 
height, and gender. 
15 Cardiac output can be further adjusted based on 

age using the formula: 

CO - 6500 ml/min *( l-0.O08*(age-30 ) ) - 
When inputting patient specific information, a choice can 
be made to further adjust the cardiac output for normal, 
20 low, and high. The cardiac output level can be estimated 
using sequential CT scanning of the aorta. 

Next, a contrast agent concentration is input and 
accepted as well as an injection method, total injection 
time, and injection rate. These values are all well 
25 known to those of ordinary skill in the art for 

particular types of CT scans. Thereafter, control is 
returned to the main program. The second step of the 
method is to call the compartment/ model subroutine. This 
subroutine, shown in Figure 8, begins with the right 
30 heart and follows the blood flow in the body model 

diagram shown in Figure 5. A standard input blood flow 
and vessel volume is assigned sequentially for each 
circle element representing a vessel compartment in 
Figure 5. Next, a standard input blood flow, a capillary 
35 volume, and an extravascular volume is assigned 

sequentially for each block element representing an organ 
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compartment in Figure 5. Thereafter^ the blood flow and 
volume of each vessel and organ compartment Is adjusted 
by the program to be proportional to the ratio of the 
cardiac output and blood volume of the patient as 
5 compared to the standard, as calculated In the 
patient/contrast subroutine. 

In the next step, a differential equation 
describing the contrast agent transport In each vessel 
compartment, as derived above, is assigned. If the 

10 element is a vessel compartment, a differential equation 
describing contrast agent transport is assigned. If the 
element is an organ compartment, two differential 
equations describing both contrast agent transport in the 
Intravascular compartment and in the extravascular 

15 compartment are assigned. Thus, each element in the 
cardiovascular system Is assigned sequentially a 
differential equation. Control is then returned to the 
main program. 

The next step in the method is to solve the 

20 differential equations which were assigned in the 

compartment /model subroutine to obtain the organ specific 
concentration. The differential equations are solved 
with numerical integration programs of the 5th-order 
Runge-Kutta method to compute contrast agent 

25 concentration as a function of time for each compartment. 

The concentration of contrast agent in an organ is 
defined as the ratio of total contrast mass at a specific 
time to the total volume of the organ. Contrast 
concentration is converted to CT enhancement in 

30 Hounsefield Units (HU) using the ratio 1 milligram I/ml = 
25 HU. The relationship between CT enhancement In HU and 
concentration of contrast agent in mg/ml depends upon 
multiple factors including the type of contrast agent, 
the surrounding tissue and other factors related to the 

35 CT scanner such as peak kllovolts used (kVp). The assumed 
relationship of 1 mg/ml equals 25 HU was arrived at 
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through an experiment comparing CT attenuation and 
contrast concentration. 

In that experiment, loversol-320 (I) was diluted 
with saline to generate various concentrations ranging 
5 from 0 to 30 mg/ml. Fifty ml deposits of the solutions 
were placed in plastic Jars and scanned with a Siemens 
Somatom Plus CT scanner using standard abdomen and chest 
settings of 120 kV^ and 137 kVp. CT attenuation was 
recorded by placing a 1.5 centimeter circular region of 
10 interest in the center of each jar on each image. 

Enhancement was computed as the difference between CT 
attenuation in each Jar and the CT attenuation in a Jar 
filled with normal saline. Figure 9 is a graph showing 
the recorded enhancement levels ranging from 8 to 800 HU 
15 for concentrations ranging from 0 to 30 mgl/ml at each of 
120 JcV and 137 kVp. When a linear relationship was 
assumed, an increase in concentration by 1 mgl/ml yielded 
an approximate increase in contrast enhancement of 25 HU. 

The last step in the method shown in Figure 6 is 
20 providing a display of the enhancement pattern of the 
vessels and organs of interest as a function of time. 
This can be through either a data stream or a graph. 

To gauge the accuracy of the invention, simulated 
graphs were generated for a hypothetical patient using 
25 different injection protocols. These simulated graphs 
were compared to empiric graphs representing actual 
enhancement level measurements in patients who had 
undergone contrast enhanced CT scans. The empiric graphs 
represent an average of the recorded enhancement levels 
30 in the aorta and liver from three groups of 25 to 28 
patients for the injection protocols listed in Table 6 
below. Each injection consisted of 125 milliliters of 
loversol-320. The data used to create the empiric 
enhancement graphs was collected in an unrelated 
35 experiment regarding enhancement levels and both 

unlphasic and btphasic injection protocols. A biphaslc 
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injection uses two injection rates during the injection 
time. A uniphasic injection uses one injection rate 
during the injection time. 

The simulated graphs represent contrast 
5 enhancement for each of the three protocols in Table 6 
based on a hypothetical patient whose weight equalled the 
average weight of the corresponding empirical group of 
patients. Thus, each point on the empiric graphs 
represents an average of a wide range of empirical 
10 enhancement values while each point in the simulated 
graphs represents a single enhancement value for a 
hypothetical patient. Figure 10 shows a simulated graph 
100 and an empiric graph 102 for the blphasic-low flow 
rate injection protocol shown in Table 6. The 
15 hypothetical patient, whose enhancement levels are 

represented in the simulated graph 100, had an assumed 
body weight of 158 pounds- This assumed body weight was 
equal to the average body weight of the 28 patients whose 
actual mean enhancement levels are represented by the 
20 empiric graph 102. 

Figure 11 shows a simulated graph 104 and an 
empiric graph 106 for the uniphasic-low flow rate 
injection protocol In Table 6. The hypothetical patient, 
whose enhancement levels are represented in the simulated 
25 graph 104, had an assumed body weight of 171 pounds. 
This assumed body weight was equal to the average body 
weight of the 25 patients whose actual mean enhancement 
levels are represented by the empiric graph 106. 

Figure 12 shows a simulated graph 108 and an 
30 empiric graph 110 for the unlphasic-high flow rate 

injection protocol In Table 6. The hypothetical patient, 
whose enhancement levels are represented in the simulated 
graph 108, had an assumed body weight of 177 pounds. 
This assumed body weight was equal to the average body 
35 weight of the 27 patients whose actual mean enhancement 
levels are represented by the empiric graph 110. 
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The simulated and empirical contrast enhancement 
graphs were compared according the maximum enhancement 
level of each graph and the percent difference between 
the graphs. The simulated graphs were In good agreement 
5 with the empiric graphs. For example. In Figure 10, for 
the blphaslc-low flow rate Injection protocol the 
simulated maximum aortic enhancement was 142.7 HU while 
the empiric maximum aortic enhancement was 163.4 HU. 
Also in Figure 10, the simulated maximum hepatic 
10 enhancement was 53.8 HU while the empiric maximum hepatic 
enhancement was 55.5 HU. 

In Figure 11, for the uniphasic-low flow rate 
injection protocol, the simulated maximum aortic 
enhancement was 220.4 HU while the empiric maximum aortic 
15 enhancement was 205.8 HU. Also in Figure 11. the 

simulated maximum hepatic enhancement was 63-8 HU while 
the empiric maximum hepatic enhancement was 59.8 HU. 

In Figure 12, for the uniphasic-high flow rate 
injection protocol the simulated maximum aortic 
20 enhancement was 321.3 HU while the empiric maximum aortic 
enhancement was 313.7 HU. Also in Figure 12, the 
simulated maximum hepatic enhancement was 63.6 HU while 
the empiric maximum hepatic enhancement was 60.8 HU. 

The total mean difference in maximum enhancement 
25 between the simulated and empiric graphs was 7.4 percent 
for aortic enhancement and 4.8 percent for the hepatic 
enhancement. As can be seen in Figures 10, 11 and 12, 
the simulated and empiric graphs were also nearly 
Identical in variation over time. Specifically, the 
30 average enhancement difference between the simulated and 
empiric graphs for all three protocols In Table 6 was 
11.6 percent for aortic enhancement and 12.7 percent for 
hepatic enhancement. 

It Is well known that body weight is one of the 
35 patient variables which most drastically affects contrast 
enhancement. To confirm the functionality of the 
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invention, the effect of body weight on contrast 
enhancement was simulated in a hypothetical patient. 
Figure 13(a) shows simulated aortic enhancement graphs 
and Figure 13(b) shows simulated hepatic enhancement 
5 graphs for uniphasic-high injection protocol in an adult 
male with a fixed height (5' 8") and body weights of 110, 
150, 200 and 250 pounds. The simulated graphs 
demonstrate that contrast enhancement was greatly 
affected by body weight. For example, in Figure 13(a), 
10 the peak aortic enhancement in a subject weight 110 

pounds was more than twice that in a subject weighing 250 
pounds. However, as expected, the timing of the aortic 
and hepatic peaks did not vary significantly because 
alteration in the cardiac output was compensated by 
15 alteration in the blood and body fluid volume. The 

simulated graphs in Figure 13 correlate well with empiric 
observations in patients showing an inverse relationship 
between body weight and contrast enhancement. 

In the patient /contrast subroutine, selection of a 
20 permeability factor between 2 and 10 is required, as 
explained, infra. However, of the variables used to 
construct the cardiovascular model, the least is known 
about the transcapillary permeability. Permeability 
varies from organ to organ and depends, in part, on the 
25 substance being transferred. Organs with discontinuous 
capillaries such as the liver, spleen and bone marrow 
have relatively high permeability. Fenestrated 
capillaries in the kidney and intestines have 
intermediate permeability. Continuous capillaries in the 
30 heart muscle and skin have smaller pores and thus lower 

permeabl 1 1 ty . 

Although some general information about 
permeability is known, knowledge about specific 
transcapillary permeability is limited. For example, the 
35 size of the contrast substance is one of the most 
important properties in determining the rate of 
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transcaplllary exchange. Permeability for different 
substances will vary according to each substance's 
molecular weight. Most nutrients and metabolites 
including glucose (mw-180) and sucrose (mw«342) are quite 
5 readily diffusible. 

When transcaplllary exchange occurs slowly 
relative to the blood flow rate, it is primarily 
diffusion- limited- Conversely, if transcaplllary 
exchange occurs rapidly relative to the blood flow rate, 
10 it is primarily flow-limited. lodinated contrast agents 
consist of relatively small molecules with molecular 
weights between 800 and 1600. Such contrast agents are 
distributed rapidly and extensively outside the blood 
vessels to the entire extracellular fluid within a few 
15 minutes of injection and are highly diffusible. 
Therefore, in the model, it was assumed that; the 
"transport of contrast agents to be mostly flow-limited 
and this assumption was applied equally to every organ in 
the cardiovascular model. 
20 Permeability (P) and transfer area (S) are usually 

treated as a unit because of the difficulty evaluating 
them separately without very detailed anatomical 
Information. The permeability- surface area product (PS) 
is referred to as the "capillary transport coefficient." 
25 The magnitude of PS in an organ is frequently expressed 
relative to the blood flow rate, Q- If PS/Q is larger 
than 1, the transport is flow-limited. If PS/Q is less 
than 1, it is diffusion-limited. In an effort to 
determine acceptable PS values in the model, simulated CT 
30 enhancement graphs were generated for several different 
PS/Q values. The simulated graphs are shown in Figure 
14(a) for aortic and Figure 14(b) for hepatic for PS/Q 
values equal to 0.1, 1, 2, 20 and Infinity. 

Simulated graphs were also generated assuming no 
35 transcaplllary barrier between the capillary and 
extracellular spaces, i.e., a single compartment 
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representing each organ. The simulated CT enhancement 
graphs generated by the Invention with PS/Q=20 closely 
approach those obtained by ignoring the transcapillary 
barrier. Thus, this PS/0 value is near the upper limit 
of flow-limited capillary transport. The simulated 
graphs shown in Figures 14(a) and 14(b), when compared 
with empiric graphs, confirm that the transport of 
contrast agent follows a flow-limited process, especially 
in richly perfused tissues. 

Figures 15 and 16 show sample aortic and hepatic 
enhancement levels generated by the invention in data 
format and graph format, respectively. Operation of the 
invention can be best understood by referring to Figure 
15. Prior to performing a CT scan, an operator inputs 
15 the patient specific information, such as height weight 
and cardiac output, and an injection protocol into the 
program in accordance with the above description of the 
invention. The program then generates output data 
showing the predicted organ specific enhancement values 
20 as a function of time. The output data can take the form 
of a data stream as shown in Figure 15 or a graph as 

shown in Figure 16. 

The operator views the data initially to determine 
whether the proposed injection protocol will result in an 

25 acceptable enhancement level for an acceptable duration. 
If the data shows that the desired enhancement level will 
never be reached, or will not be sustained for a 
sufficient length of time, the operator chooses a 
different injection protocol and then reruns the program 

30 until a satisfactory predicted enhancement level is 
obtained . 

After the operator obtains an output showing an 
acceptable predicted enhancement level and duration, the 
operator then selects a scan start time and duration, 
35 including an appropriate collimatlon thickness and table 
speed. in the alternative, all or a portion of the 
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selection can be performed by the computer. This 
Information is then input into the CT scanner, if 
obtained off-line from the CT control computer, and the 
scan is then executed. For example, assuming a threshold 
5 hepatic enhancement level of 50, the data of Figure 15 
shows that the threshold enhancement level is not reached 
until .64 minutes after the injection of contrast agent 
into the patient. In addition, the data shows that the 
threshold enhancement level will be maintained for 
10 approximately 1.7 minutes. Using this information, an 
operator inputs the scan start time, scan duration, 
coll imat ion thickness and table speed into the CT scanner 
and thereafter performs the scan on the patient. In the 
alternative, computer software can be implemented to 
15 automatically transmit the output information directly 
into the CT scanner. 

The program of the invention gives an output which 
includes a predicted enhancement level in the tissue of 
interest as a function of an elapsed time after 
20 injection. As discussed above, the enhancement threshold 
is a level of tissue enhancement below which results ±n a 
poor quality scan. The scan duration is the time between 
starting the scan and ending the scan. In order to 
optimize the scan, the tissue scanned must maintain an 
25 enhancement level equal to or greater than the threshold 
enhancement level for the entire scan duration. 

The inventors herein have further expanded the 
invention by providing a means for analyzing the 
generated predicted enhancement function (enhancement 
30 level with respect to time) to determine if the predicted 
enhancement level in the tissue to be scanned 
sufficiently meets the criteria required for an optimum 
scan and providing a means for adjusting the injection 
protocol until the output is acceptable. 
35 In the preferred embodiment, preferably 

Implemented in a computer program, if the output 
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predicted enhancement: function indicates that the tissue 
enhancement level will never attain the desired 
enhancement threshold or that the tissue enhancement 
level will not be maintained above the desired threshold 
5 for a period of time equal to or exceeding the scan 

duration, the program provides a means for adjusting the 
injection protocol. 

Referring to Figure 18, the output of the 
preferred embodiment is a graph or curve 150 showing a 

10 predicted enhancement level as a function of time 

superimposed on a line 152 representing the enhancement 
threshold. When the peak enhancement, the highest level 
reached by curve 150, is below the desired threshold, or 
when the time interval (B-A) is shorter than the desired 

15 scan duration, enhancement level must be raised to obtain 
an acceptable scan* In the graph of Figure 18, 
enhancement curve 150 exceeds the enhancement threshold 
for time period (B - A), 

If the peak enhancement does not reach the 

20 threshold or if the time interval (B-A) is not equal to 
or greater than the desired scan duration, the preferred 
embodiment provides the operator two options to increase 
the predicted enhancement level. It is known in the art 
that increasing contrast volume, flow rate, or 

25 concentration increases the level of enhancement. 

Because enhancement level is a function of the amount of 
contrast transported through a particular tissue over a 
given time period, increasing flow, volume and 
concentration all result in Increased enhancement levels. 

30 In the majority of the injector systems, volume and flow 
rate are easily adjusted through the injector system 
controls. In addition, most medical facilities have a 
limited number of differing concentrations. Therefore, 
the most practical adjustments to adjust enhancement 

35 levels are to volume and flow rate. However, adjusting 
concentration is also acceptable to raise enhancement. 
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The steps of this aspect of the preferred 
embodiment are shown In the flow chart of Figure 21. 
After receiving an Input of a predicted contrast 
enhancement function, an enhancement threshold and a scan 
5 duration, the program determines if the predicted 

enhancement function exceeds the threshold for a time 
period greater than the scan duration. If not, the 
operator has the option of either (i) increasing the 
contrast volume or (11) increasing the flow rate. After 
10 receiving new volume, flow rate or both, the program 

determines if the maximum allowable flow rate or maximum 
allowable volume have been reached. It is known in the 
art that Injection flow rates and volumes have maximum 
limits above which safety concerns are implicated. These 
15 limits vary depending on many factors including the 

particular patient, the contrast agent and the particular 
procedure involved. Therefore, the present invention 
provides for the input of maximum values for flow rate 
and volume. When either the maximum allowable contrast 
20 volume or injection rate is reached without achieving the 
threshold enhancement level, the program continues to 
allow increasing the other until both maximums are 
reached. After input of new values, the program iterates 
the steps of updating the predicted enhancement function 
25 based on new values. The contrast volume or the flow 

rate is progressively increased in this fashion until the 
time interval (B-A) becomes equal to or greater than the 
requested scan duration. 

in the preferred embodiment, once the volume and 
30 flow rate have reached the maximum values, and the 
predicted enhancement function does not exceed the 
threshold for a time period greater than or equal to scan 
duration, the program notifies the operator that a new 
enhancement threshold must be chosen. 
35 Although the preferred embodiment allows the 

operator to input different flow rates and volumes, the 
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enX:±T& process could be automated and performed by a 
computer processor. For example, as an alternative to 
the operator selecting the flow rate and volume to obtain 
an acceptable enhancement level, a linear bisection 
5 method, or other known mathematical process, could be 
programmed into the computer to solve for a convergence 
point by reducing the differences between and updating 
two boundary values. 

Although a variance in cardiac output may affect 
10 the level of enhancement in a given tissue, the inventors 
herein have discovered that a change in cardiac output 
will more dramatically affect the time at which a 
particular level of enhancement will be achieved. 
However, the program selects a cardiac output believed to 
15 be most closely associated with the patient and uses that 
value In its computation of the predicted enhancement 
function. The program also provides for operator Input 
of alternative values for cardiac output with each 
resulting enhancement function being tested using the 
20 above described method. In this way, an operator can be 
certain that a specific injection protocol will result in 
a predicted enhancement function which exceeds the 
desired threshold for the entire scan duration regardless 
of the cardiac output of the patient. 
25 As shown in the flow chart of Figure 21, when the 

predicted duration of the enhancement (B-A) is greater 
than the requested scan duration, the method of the 
present invention allows the operator to select from the 
options of (1) reducing the volume of the contrast medium 
30 or (ii) reducing the injection rate; or (iii) maintaining 
the current enhancement level and searching for the 
optimal scan interval within A and B. Options (1) and 
(11) may be used for a more efficient scan or for reasons 
related to a patient's medical history. For example, by 
35 reducing the volume of the contrast medium, costs are 
saved and the patient need not be given unnecessary 
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additional contrast to create an acceptable scan. 
Because the contrast medium might have side effects on 
the patient, reducing the amount of contrast may be 
desired to limit the amount of contrast a patient miast 
5 receive in order to undergo a successful CT scan. 

AS shown in the flow chart of Figure 21, the 
process steps of decreasing the voUane and/or rate are 
virtually identical to the method of increasing these 
values, except for the direction of adjustment, and the 
10 discussion above related thereto is equally applicable 
here. It is also foreseeable that the process steps for 
increasing and decreasing the flow rate and/or volume 
could both be utilized In one scan procedure, if for 
example, an adjustment In one direction resulted in too 
15 great a change in enhancement level or, after lowering 
the threshold, the enhancement level is predicted to 
exceed the revised threshold for an excessive period. 

In the alternative, or If the above concerns are 
not a factor, as mentioned above in option (lil) the 
20 preferred embodiment allows the option of maintaining the 
selected injection protocol and determining the optimum 
interval between time A and time B which is equal to the 
scan duration and during which the tissue enhancement 
level is the greatest. The method described initially 
25 herein was a significant improvement over the prior art 
m that It allowed prediction, prior to injection, of 
whether the enhancement level would exceed the threshold 
value or whether modifications were required In the 
various input parameters to obtain a tissue enhancement 
30 level which exceeded the threshold value for the scan 

duration. In this fashion, an operator can determine the 
proper delay after initiating the injection to begin 
scanning as well as an optimum scan duration. This 
information also enabled the operator to change various 
35 scan parameters, including table speed and collimatton 

thickness in order to achieve a scan during the time when 
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the tissue enhancement level exceeded the threshold 
level. 

Under certain circumstances, the predicted tissue 
enhancement function exceeds the threshold enhancement 
level for a time period greater than the scan duration. 
The operator could choose arbitrarily at what point after 
the tissue enhancement level exceeded the threshold to 
begin the scan, as long as the scan would be completed 
before the tissue enhancement level decreased below the 
threshold level and obtain an acceptable scan. Building 
further on the invention, the inventors herein have 
provided a means for selecting a scan start time to cause 
the scan to take place during the optimum temporal window 
between time B and time A if the duration of predicted 
enhancement above threshold exceeds the scan duration by 
more than a predetermined amount of time. As shown in the 
flow chart of Figure 21, when the time interval (B-A) is 
equal to (or approximately 10 % greater than) the 
expected scan duration, point A becomes the onset of 
scanning. Thus, the optimal temporal window of scanning 
begins at time A and ends at time B. If, however, the 
time interval (B-A) is significantly greater than the 
expected scan duration plus 10%, the invention selects an 
optimal temporal window within B-A by maximizing the 
predicted enhancement available. This is best understood 
by referring to Figure 19. 

in Figure 19, the enhancement curve or function 
150 of the tissue to be scanned is displayed along with a 
desired threshold of enhancement 152 on a graph of 
30 enhancement versus elapsed time after injection. As 

shown, the enhancement level of the tissue to be scanned 
increases from level 0 at time 0 to a level equal to the 
desired enhancement threshold at time A. The tissue 
enhancement level continues to rise to a peak enhancement 
35 level above the threshold and then decreases to again 

equal the desired threshold at time B. After time B, the 
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tissue enhancement: level continues "to decrease below the 
threshold. As Is known in the art, a unlphaslc Injection 
results in a single peak enhancement level being attained 
as shown in Figure 19. Thus, the duration of predicted 
5 enhancement above threshold is represented as (B-A). 

Referring still to Figure 19, the method of the 
present invention selects two points (C and D) whose 
difference (D-C) is equal to the scan duration and 
between which the area under the enhancement curve above 

10 the threshold (AUG) is maximized. As executed in one of 
the steps of the flow chart of Figure 21, a set of AUC's 
are calculated as point C Incrementally advances with a 
fixed interval (D-C) to determine the optimum scan 
interval- The first AUC is calculated when C coincides 

15 with A. This would correspond to beginning the scan at 
the time the tissue enhancement level first equals the 
threshold value. The last AUC Is calculated when D 
coincides with B. This would correspond to ending the 
scan when the tissue enhancement level equals the 

20 threshold level for the second time. After calculating 
the set of AUC's, choosing the maximum AUC sets the 
optimal scan start time (C) and scan end time (D). 

The number of AUC's in the set is dependant on the 
Increment with which C Is advanced. Of course, the 

25 smaller the increment of advancement, the more accurate 
on average the time Interval predicted for optimal 
scanning will result. The Inventors herein have found 
that calculating the AUC's while incrementally advancing 
C in units of 1 second Is satisfactory. When the 

30 enhancement curve is unlphaslc^ a plot of the set of AUCs 
demonstrates a similar distribution and the optimal 
scanning Interval contains the peak enhancement point. 
Thus, in a unlphasic scan the maximum AUC should be known 
once the calculated AUC begins to decrease from a peak 
35 value. 
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Unlike the prior art, the Invention inltlall 
disclosed herein allows the operator to predict, prior to 
starting an injection, whether tissue enhancement will 
successfully attain a threshold, whether the tissue 
5 enhancement level will be maintained above the threshold 
for the entire scan duration, and the proper scan delay 
to allow the scan to begin at a time when the tissue 
enhancement level exceeds the threshold value. The above 
described improvements enhance the invention by providing 
10 a means for optimizing the injection protocol to obtain 
acceptable enhancement levels and optimizing the scan 
start time when the period of acceptable enhancement is 
predicted to be greater than the scan duration. 

These methods and those initially described are 
15 significant improvements over the prior art and allow 
prediction well within acceptable limits of accuracy. 
However, the inventors herein have further improved the 
invention by providing for even more accurate 
predictions. As explained above, the tissue enhancement 
20 level is directly related to the volume and concentration 
of contrast in the tissue to be scanned when the scan 
takes place. Because the contrast agent is distributed 
throughout the patient by the cardiovascular system, the 
amount of contrast in a given tissue at a given time is 
25 related to various patient specific parameters which 

affect contrast transport throughout the patient. These 
include height, weight, gender, age, and cardiac output. 

Cardiac output, unlike height, weight, gender and 
age is not readily measured. As explained above, many 
30 factors including disease status or prior heart failure 
may affect cardiac output. In addition, a patient does 
not maintain the same cardiac output during his or her 
entire lifetime. 

inputting the correct cardiac output into the 
35 model is necessary for accurate prediction of optimum 

scan delay. For example, if a patient has a poor cardiac 
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output and the operator uses a standard cardiac output 
because the cardiac status ±s unknown, the predicted 
optimum time to begin scanning a tissue will not coincide 
with the actual optimum time to begin scanning the 
5 tissue. Although the predicted threshold level will 

eventually be met, the maximum period of enhancement may 
still not coincide with the period of scanning due to the 
input of an Inaccurate cardiac output. 

To account for varying cardiac output, the 
10 intention initially described provided for the input of 
different cardiac output values Into the mathematical 
model of the cardiovascular system. m addition, the 
invention Initially described provided for the operator 
to choose several alternative values for cardiac output 
15 and generate a family of predicted tissue enhancement 
functions or curves corresponding to the alternative 
values- By predicting and analyzing the family of curves 
prior to injection, the operator could be certain that 
the tissue enhancement level would meet or exceed the 
20 threshold enhancement regardless of the cardiac output. 

In the later described invention, each family 
member has a slightly different cardiac output status so 
that the entire family represents the predicted tissue 
enhancement functions for the entire spectrum of cardiac 
25 output status. Increasing the number of family members 
Increases the accuracy of the prediction for a given 
patient. For example varying standard cardiac output by 
10% between each family member from 10% to 110% will give 
the predicted tissue enhancement level in ten patients 
30 having Identical body habitus, except for cardiac output, 
and the predicted enhancement functions will reflect the 
difference caused solely by the 10% difference in cardiac 
output between each patient. 

The Invention Initially described provides for 
35 taking actual measurements of enhancement after 

initiating the injection and comparing the family of 
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enhancsement curves to the actual enhancement curve to 
determine which family member most closely resembles the 
actual enhancement levels. In that way, it provides a 
means for determining early in the scan, before the 
5 threshold level had been reached, whether the scan 

parameters were appropriate and allows for adjustment if 
necessary . 

Although this is a significant improvement over 
the prior art, the comparison of the family of curves to 

10 the actual enhancement measurements in the tissue to be 
scanned must be performed quickly to allow time to adjust 
the scan parameters if necessary. Therefore, the later 
described invention uses a combination of a predicted 
enhancement function, in a region of interest, such as 

15 the aorta, prior to injection, and measurements of actual 
enhancement levels, after initiation of injection, to 
calculate a correction factor (such as proper cardiac 
output status) to be used by the model when predicting a 
tissue enhancement function for the tissue to be scanned, 

20 such as the liver. The preferred embodiment provides for 
sequential low-dose pre-scanning of the aorta after 
injection to calibrate the mathematical model to unknown 
or difficult to measure specific patient parameters such 
as cardiac output. Because the contrast reaches the 

25 aorta quicker than the liver, using enhancement level 

measurements in the aorta as feed back increases the time 
allowed to make corrections to the scan parameters before 
the required onset of hepatic scanning. 

The preferred embodiment is implemented in a 

30 computer program and can be implemented in a stand alone 
computer, a computer included in a CT scan machine or a 
computer included in an Injector system. Moreover, the 
present Invention could be implemented in a CT scan 
system which includes an injector and a CT scan machine 

35 both controlled by the same computer, in the preferred 
embodiment, the computer program contains a mathematical 
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model of the patient's cardiovascular system. The 
details of the mathematical model are fully explained 
above . 

The computer program accepts input values for 
5 those parameters in the patient and in the injection 
protocol which affect contrast transport through the 
cardiovascular system. These include patient age, 
gender, height, weight, cardiac output and injection flow 
rate, volume, concentration, phase and scan duration. 
10 The program accepts the inputs and generates a predicted 
enhancement level as a function of an elapsed time after 
injection for both aortic and hepatic enhancement. 
Figure 20 shows a table giving the data output from the 
program for a particular patient having a particular body 
15 habitus and presuming a standard cardiac output. The 
operator has already adjusted the injection protocol 
using the method set forth above to ensure that the 
predicted enhancement function will exceed the threshold 
for a length of time exceeding the scan duration. The 
20 columns of information displayed in Figure 20 are as 

follows: the left most column headed "time" represents 
elapsed time from the start of injection; the next column 
to the right displays the calculated predicted hepatic 
enhancement level in Hounsefield Units; the next column 
25 is the difference between the predicted enhancement level 
and the preselected threshold chosen as 50 Hounsefield 
Units; and the next column displays the calculated AUC, 
with each entry being aligned with its corresponding scan 
start time. For example, the first entry of 9.6 
30 corresponds to a scan start time of 40 seconds and a scan 
end time of 70 seconds. Using the methodology of the 
present invention, the maximum AUC is readily identified 
as 263.8 which corresponds to a scan start time of 60 and 
a scan end time of 90. The start scan and end scan times 
35 for the optimal temporal window are identified in the 
last column of the table of Figure 20. 
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TO enhance the prediction of scan start time, the 
method of the present invention can be used to update or 
confirm the output of the table in Figure 20. First, a 
base line scan is performed at a region of interest prior 
5 to the initiation of the injection so as to enable a 
calculation of the actual enhancement level for the 
region of Interest. The base line scan, as is Known In 
the art, is a low-dose or partial scan in which the x-ray 
dose is reduced sxibstantially compared to a typical scan 
10 and views may be acquired through less than a full 
revolution of the gantry. The x-ray dose is thus 
considerably less than a normal image scan, but 
nevertheless, a slice image may be reconstructed. 

In the preferred embodiment, an actual aortic 
15 enhancement function Is compared with predicted aortic 
enhancement functions, generated using different cardiac 
outputs, to calculate a correction factor to be applied 
to the model before predicting a hepatic enhancement 
function. After the correction factor is applied, the 
20 predicted hepatic enhancement function generated is more 
accurate than the predicted hepatic enhancement function 
generated without the correction factor. Although the 
region of interest monitored can be the tissue to be 
scanned, it is preferable to monitor a region of interest 
25 distinct from, and one which will provide a measured 

response faster than, the tissue to be scanned to allow 
for maximum time to calibrate the mathematical model to 
the particular patient and update the scan parameters. 
In the preferred embodiment, a region of interest is 
30 selected that can be monitored for actual enhancement and 
analyzed sufficiently before the tissue to be scanned 
attains a threshold enhancement level. This allows the 
program to use the feedback from the monitoring to 
enhance the accuracy of the model which is then used to 
35 predict a tissue enhancement function. The more accurate 
tissue enhancement function can then be used to select 
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optimum scan parameters well be£ore "the onset of 
scanning • 

After the operator has selected an appropriate 
region of Interest and performed a base line scan of the 
5 region of Interest^ the Injection Is Initiated according 
to the selected Injection protocol. After the injection 
Is initiated^ and as the injection is being administered, 
the CT scan machine is used to monitor the region of 
interest enhancement level using low-dose pre-scan 
10 x-rays. The information from the pre-scan monitoring is 
used by the program to generate an actual aortic 
enhancement function as shown in Figure 23. As In Figure 
23, the predicted 160 and actual 162 regional enhancement 
function can be displayed on the same graph for 
15 comparison. The data generated by the pre-scan 

monitoring can also be used to display the contrast 
enhancement level in a chart or to reconstruct an actual 
image of the region of interest being monitored. 

After sufficient time has elapsed, the predicted 
20 regional enhancement function can be compared to the 
actual regional enhancement function generated by the 
low-dose pre-scannlng of the region of interest. The 
results of the comparison can be used to calibrate the 
mathematical model for factors such as cardiac output 
25 before generating the predicted tissue enhancement level 
for the tissue to be scanned. Referring to Figure 23, 
the predicted aortic enhancement function 160 is shown 
along with the actual aortic enhancement function 162 
from pre-scan monitoring. 
30 Also shown in Figure 23 are the originally 

predicted hepatic enhancement function 164 and an updated 
or calibrated predicted hepatic enhancement function 168 
corresponding to the two aortic curves. As can be seen, 
due to a measured decreased cardiac output in the 
35 patient, the onset of hepatic scanning must be delayed 
longer than originally predicted. It is this feedback 
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that: ±s used by the present invention to calibrate or 
fine tune predicted hepatic scanning based on the actual 
measurement of aortic enhancement levels with low-dose 
pre- scanning . 

5 The slope of the measured aortic enhancement curve 

162 may be calculated at a predetermined time and 
compared to the calculated scope of the predicted aortic 
enhancement curve 160 to determine the difference between 
the predicted and actual patient aortic output. However, 

10 this slope comparison has been unreliable as it has been 
observed that early aortic enhancement measurement is 
frequently pulsatile and noisy. The inventors herein 
have discovered that by graphing the actual enhancement 
as a function of time elapsed after injection (the 

15 enhancement curve 162) and measuring the area under the 
enhancement curve (AHC) after a predetermined time 
interval, this calculated AHC provides a much more 
reliable indicator of the patient's aortic output. 

The method of the preferred embodiment plots the 

20 actual measurements of aortic enhancement to represent 
the actual enhancement level as a function of time and 
calculates the area under the curve (AHC) of the actual 
aortic enhancement function at a predetermined time. 
This AHC can compared to the area under the curve of the 

25 predicted regional enhancement function to gauge the 

accuracy of the model and calculate a correction factor. 
The present invention is very useful in determining 
predicted enhancement levels in specific tissue of a 
patient when not all of the patient specific parameters 
30 which affect transport are known. For example, as 

discussed above, cardiac output is a patient specific 
parameter which is not readily measurable. However, the 
differences between a predicted aortic enhancement 
function and an actual measured enhancement function may 
35 be analyzed to determine the cardiac output of a patient. 
It is known that the delay between aortic and hepatic 
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enhancement represents the time required to distribute 
contrast medium from the aorta to the liver and Is 
proportional to the cardiac output. The slower the 
cardiac output, the longer the delay between the time 
5 that contrast medium Is delivered to the aorta and the 
time It Is delivered to the liver. Thus, with less than 
standard cardiac output the onset of hepatic scanning 
must also be delayed so as to coincide scanning with peak 
enhancement . 

10 The method of the present invention was used to 

exhibit the effect of cardiac output on contrast 
enhancement for a hypothetical adult male with a fixed 
height (5 ft. 8 In.) and body weight (150 lbs) subjected 
to uniphaslc-hlgh injection protocol. The cardiac output 
15 specified for the model was varied by multiplying the 
standard cardiac output by 0.25, 0.50, and 2.0. Four 
enhancement curves were generated for predicted aortic 
enhancement as shown in Figure 22(a) and four enhancement 
curves for predicted hepatic enhancement as shown in 
20 Figure 22(b). As can be seen from Figure 22(a), as the 
cardiac output decreases, the time delay to the peak 
enhancement Increases in both aortic and hepatic 
enhancement curves. The peak aortic value Increases with 
reduced cardiac output, while in Figure 22(b) the plateau 
25 of peak hepatic enhancement is prolonged. 

The method of the present invention is practiced 
es follows. A look up table can be constructed using the 
computer program of the present invention to generate 
different outputs of predicted regional enhancement 
30 levels In a region of Interest based upon different 

cardiac outputs. An example of such a table is shown in 
Figure 24. As shown in the table, predicted aortic 
enhancement levels are calculated at 5 second intervals 
for high cardiac output (200%), standard cardiac output 
35 (100%), reduced cardiac output (75%) and low cardiac 
output (50%). The area under the aortic curve is also 
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calculat:ed at intervals of 5 seconds for each different 
cardiac output. The last coltimn of Figure 24 shows an 
example of a list of actual aortic enhancement levels 
measured from low dose pre- scanning and areas under *the 

5 actual regional enhancement curve for those enhancements. 

A predetermined time after Injection Is chosen, 
for example 20 seconds, shown In the last column of the 
table In Figure 24, the actual AHC at time equal 20 
seconds is calculated. The operator (or the computer) 

0 can then compare the actual AHC at 20 seconds with the 
predicted AHCs in the first four columns to determine 
which column most accurately predicts an AHC of 776 
HU*sec. at time 20 seconds. As can be seen from the 
table, the column having cardiac output equal to 75% most 

5 closely matches the AHC at 20 seconds (AHC»768 HU*sec). 
Once the cardiac output is determined in this fashion, 
the program can calculate a predicted hepatic enhancement 
function using cardiac output of 75%. This will allow 
the prediction of hepatic enhancement with much more 

O accuracy because the model is more closely calibrated to 
the patient specific parameters which affect contrast 
enhancement • 

The table of Figure 20 was previously explained to 
determine the optimxim onset of hepatic scanning for a 

5 patient with a specific body habitus and assuming a 

standard cardiac output. Using the results of the aortic 
monitoring, shown in the table of Figure 24, a revised 
table was generated for the same patient. As described 
above, the results of the table in Figure 24 reveal that 

0 a more accurate cardiac output for the patient is 75% of 
the standard. Therefore, a new table, shown in Figure 
25, was generated for predicting the optimum onset of 
hepatic scanning. As can be seen comparing Figure 20 to 
Figure 25, the optimum onset to hepatic scanning is 

5 changed from 60 seconds, which was predicted using 
standard cardiac output, to a delay of 80 seconds. 
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predlct:ed using cardiac output of 75% standard. The 
optimum scan Interval Is determined for the updated 
hepatic scanning parameters using the present invention 
by the greatest AHC for the Interval of scan duration as 
5 discussed above. Thus, the hepatic enhancement levels 
can be accurately predicted to allow the entire scan 
duration to take place during an interval of maximum 
enhancement even if cardiac output cannot be readily 
measured . 

10 Furthermore, the operator, or computer, has 

sufficient time to modify the scan parameters, such as 
scan delay, because the calculation of cardiac output was 
performed at 20 seconds after injection initiation and 
the onset of hepatic scanning should not occur until 

15 approximately 80 seconds after initiation of injection. 
Although a greater delay before calculating the area 
under the actual aortic enhancement curve could provide a 
more accurate outcome, the inventors have found that 
measuring aortic output at 15-20 seconds after initiation 

20 is acceptable and allows sufficient time to update the 
scan parameters. 

An alternative embodiment of the method of the 
present invention is shown in the flow chart of Figure 
26. This embodiment, preferably embodied in a computer 

25 program, calls for the creation of an interpolating look- 
up table showing the times (tl) that different cardiac 
outputs reach a predetermined AHC and showing predicted 
optimum hepatic scanning times (tl*) correlated to each 
(tl). The input to this embodiment is the actual level 

30 of aortic enhancement which is being attained in the 
patient. 

The program uses the input to graph the actual 
enhancement function and compute at regular intervals the 
area under the actual enhancement curve (AHC). When the 
35 computed actual AHC is equal to the AHC from which the 

table was generated, the time expired after the injection 
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Initiation Is recorded. The recorded time Is then used 
vrlth the table to find the tl closest to the recorded 
tine. Once the corresponding tl Is found In the table, 
the corresponding optimum time for onset of hepatic 
5 scanning (tl*) Is given. Using the table, the operator 
can thus obtain the predicted delay prior to beginning 
the onset of hepatic enhancement well before the time to 
begin the hepatic scan occurred. 

There are various changes and modifications which 
10 may be made to the Invention as would be apparent to 
those skilled In the art. However, these changes or 
modifications are Included In the teaching of the 
disclosure, and It Is Intended that the Invention be 
limited only by the scope of the claims appended hereto. 
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EXHIBIT A 

TABLE 1. ESTIMATED DISTRIBUTION OF BLOOD IN VASCULAR SYSTEM OF AN ADULT HUMAN 

VOLUME 

REGION ^ I 

HEART (DIASTOLE) 360 7.2 

PULMONARY 440 B.6 
ARTERIES 130 2.6 

CAPILLARIES ISO 3.0 

VEINS 160 3.2 

SYSTEMIC 4.200 81 
AORTA ANO LARGE ARTERIES 300 6.0 

SMAU ARTERIES 100 8.0 

CAPILLARIES 300 6.0 

SMALL VEINS 2.300 16.0 

LARGE VEINS 900 18.0 

TOTAL 5^00 100 

MODIFIED FROM THE REFERENCE (NILNOR) 



TABLE 2. ESTIMATED DISTRIBUTION OF CARDIAC OUTPUT IN AN ADULT HUMAN. 

THE LIVER RECEIVES DUAL BLOOD SUPPLIES. HEPATIC ARTERY AW) PORTAL SYST EM 

BLOOD FLOW 

REGION ■L/NIM. X ~ 



UPPER EHREHITIES 325 S.O 

HEAD 975 IS.O 

CORONARY 260 4.0 

BRONCHIAL 130 2.0 

KIDNEYS 1.430 22.0 

LIVER 1.885 29.0 
HEPATIC ARTERY 455 7.0 
PORTAL 1.430 22.0 



SPLEEN/STOMACH 



430 6.6 



PANCREAS/INTESTINE . „ „15.4 

TRUNK/LOUER EXTREMITIES ^-^95 23-0 

TOTAL 100 

MODIFIED FROM REFERENCE (UAOE) 



TABLE 3. ESTIMATED DISTRIBUTION OF BODY FLUID IN AN ADULT HUNAN. 

COMPARTMENT VOLUME (LITER) 

ICF (EXCEPT RBC) 19.1 
ECF (EXCEPT PLASMA) 15.9 
BLOOD (PLASMA«RBC) S.O 



TOTAL BODY WATER ^0 



ESTIMATION BASED ON THE VOLUME OF DISTRIBUTION OF lOHEXOL (OLSSON). 
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TABLE i. ESTIMATED DISTRIBUTION OF Z BLOOD FLOU RATE AND CAPILLARY VOLUNE IN AN 
ADULT HUHAN. REGIONAL CAPILLARY VOLUHE IS CALCULATED PROPORTIONAL TO A REGIONAL 
BL€OD FLOU OF THE TOTAL 122Z f THE PORTAL SYSTEH CONTRIBUTES AN ADDITIONAL 22Z) 

~" BLOOD FLOW CAPILLARY VOLUH E 
REGION X 

UPPER EXTRENITIES S.O 12 

HEAD IS.O 37 

CORONARY 1.0 10 

BRONCHIAL 2.0 5 

KIDNEYS 22.0 54 

LIVER 29.0 71 

SPLEEN/STONACH 6-5 If 

PANCflEAS/INTESTINE J^-^ 38 

THUNKAOWER EXTRENITIES 23 0 S7 

TOTAL 100 C 122 ) 300 



TABLE S. ESTIHATEO WEIGHT. WATER CONTENT. AND FLUID VOLUHE OF VISCERAL ORGANS IN A 
70 Kg. AOULT. THE LUNG CONSISTS OF SOZ PARENCHYHA AND SOZ NON-PARENCHYHA TISSUES 
WHOSE CAPILLARY VOLUHES ARE ISO bL AND S ri.. RESPECTIvav. 



ORGAN 


WEIGHT (g) 


Z WATER 


FLUID ( iL ) 


FLUID-CAP. 


BRAIN 


l.iSO 


76 


1.102 


1.065 


HEART 


300 


79 


237 


227 


LUNG 


500*500 


79 


790 


635 


KIDNEYS 


300 


83 


219 


195 


LIVER 


1.800 


68 


1.221 


1.1S3 


SPLEEN/STOHACH 


170/150 


70 


221 


208 


PANCREAS/INTESTINE 


60/1.770 


70 


1.281 


1,213 


TOTAL 


7.000 




5.107 


1.726 



NUIFIED FRQN REFERENCES (ICRP. HAPLESON) 



TABLE 6. TESTED INJECTION PROTOCXS 



FIRST 

PROTOCOL 


FIRST 
RATE 
VOLUME 
(■L) 


SECOND 
RATE 
(ri./SEC) 


SECOND 
RATE 

VOLUME 
CIL) 


INJECTION 
TIME 
CSEC) 


NUMBER 

OF 

PATIENTS 


HEAN 

(fUNGE) 
OF WEIGHT 
(lb>) 


BIPHASIC-LOW 2.5 


50 


1 


75 


95 


28 


158 ( 100-205 ) 


UNIPHASIC-LOU 2.5 


125 






50 


25 


171 (108-211) 


UNIPHASIC-HI6H 5.0 


125 






25 


27 


177 (98-300) 
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What Is claimed Is ; 

1. A method of scanning a tissue In a patient 
using computed tomography, said patient having a 
plurality of patient specific parameters, wherein the 
tissue to be scanned Is enhanced with an Intravascularly 

5 Injected contrast agent, comprising the steps of: 

selecting an Injection protocol for the contrast 
agent ; 

predicting prior to implementing said injection 
protocol a tissue enhancement level as a function of an 
10 elapsed time after Injection based on the injection 
protocol and said patient specific parameters; 

determining a set of parameters for an optimum scan 
based on the predicted tissue enhancement level; and 
performing the scan in accordance with the pre- 
15 determined set of optimum scan parameters. 

2. The method of Claim 1 wherein the step of 
predicting the tissue enhancement level comprises the 
steps of: 

providing a mathematical model which describes 
5 contrast agent transport throughout the patient ' s 
cardiovascular system; 

Inputting into the mathematical model as said 
patient specific parameters patient specific information 
and contrast agent specific Information which affects 
10 physiological parameters of contrast enhancement; and 

using the mathematical model to compute a contrast 
agent concentration as a function of time. 

3. The method of Claim 2 wherein the step of 
providing a mathematical model Includes the steps of : 

providing a compartmental model of the human 
cardiovascular system Including vessel compartments 
5 representing blood vessels and organ compartments 
representing organs; 
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assigning at least: one differential equation 
describing contrast agent transport to each vessel 
compartment ; and 
10 assigning at least one differential equation 

describing contrast agent transport to each organ 
compartment • 

4 • The method of Claim 3 further comprising the 
steps of: 

assigning a standard input blood flow and a standard 
vessel volume for each vessel compartment; 
5 assigning a standard input blood flow, a standard 

capillary volume, and a standard extravascular volume for 
each organ compartment ; 

adjusting the standard input blood flow of each 
vessel compartment and each organ compartment as being 
10 proportional to a ratio of cardiac output of the patient 
to a standard cardiac output; and 

adjusting the standard capillary volume and the 
standard extravascular volume of each organ compartment 
as being proportional to the ratio of blood volume of the 
15 patient to a standard blood volume. 

5. The method of Claim 4 wherein each organ 
compartment is further subdivided into a capillary 
compartment and an extravascular compartment and the step 
of assigning at least one differential equation to each 

5 organ compartment includes assigning at least one 

differential equation to each capillary compartment and 
each extravascular compartment* 

6. The method of Claim 2 wherein the inputting 
step includes the steps of: 

inputting a permeability factor; 
inputting a patient habitus; 
5 inputting a patient blood volume; 

inputting a patient cardiac output; 
inputting a contrast concentration value; 
inputting an injection method; 
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inputting an injection time; and 
10 Inputting an injection rate. 

I . The method o£ Claim 6 wherein the step of 
inputting a patient habitus includes inputting the 
patient's age, gender, weight, and height. 

8. The method of Claim 7 wherein the step of 
inputting the patient blood volume includes the step of 
computing the patient's blood volume. 

9* The method of Claim 8 wherein the step of 
inputting the patient cardiac output includes the steps 
of calculating a standard cardiac output and adjusting 
the standard cardiac output proportional to a ratio of a 
5 cardiac output of the patient to a standard cardiac 
output . 

10. The method of Claim 2 further comprising the 
steps of: 

converting the computed contrast agent concentration 
into the tissue enhancement level as expressed in 
5 computed tomography enhancement units; and 

displaying the tissue enhancement level in computed 
tomography enhancement units as a function of time. 

II. A method for predicting a tissue enhancement 
level as a function of time for a specific tissue in a 
patient receiving an intravascularly injected contrast 
agent via a specified injection protocol, said method 

5 comprising the steps of: 

providing a mathematical model of the cardiovascular 
system of the patient, said model mathematically 
describing the transport of said contrast agent through 
said cardiovascular system, 
10 inputting into said model a plurality of patient 

specific parameters which impact the transport of said 
contrast agent, 

using said model to calculate a predicted 
concentration of said contrast agent in said tissue as a 
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15 function of time ±n response to said specified injection 
protocol, and 

displaying said predicted contrast agent 
concentration as a function of time as representative of 
said tissue enhancement level as a function of time. 

12. The method of Claim 11 further comprising the 
step of inputting a set of parameters into said model 
corresponding to said specified injection protocol^ said 
predicted concentrations being responsive thereto. 

13. The method of Claim 12 further comprising the 
step of comparing a mlnimiim threshold value with said 
tissue enhancement level to thereby determine a start 
time and time duration of acceptable tissue enhancement 

5 level during which a CT scan may be taken. 

14. The method of Claim 13 wherein said 
mathematical model comprises a plurality of compartments 
representative of significant organs and vessels in said 
cardiovascular system, each of said compartments being 

5 characterized by at least one differential equation. 

15. The method of Claim 14 wherein said 
mathematical model is embodied In a computer program. 

16. A computer being programmed for predicting a 
structure specific CT enhancement level in a patient for 
a given patient habitus and a specific contrast injection 
protocol comprising: 

5 a computer having a memory; 

a computer program in said memory, said program 
having means for accepting patient specific information 
and contrast specific information which affect 
physiological parameters of contrast enhancement and 
10 means for computing and outputtlng operator selected 
organ specific contrast concentration as a function of 
time. 

17. The computer of Claim 16 wherein the means for 
computing organ specific contrast concentration comprises 
a mathematical compartmental model of the human 
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cardiovascular system including vessel compartments 
5 representing blood vessels and organ compartments 
representing organs; 

at least . one differential equation describing 
contrast transport to each vessel compartment; and 
at least one differential equation describing 
10 contrast transport to each organ compartment. 

18. The computer of Claim 17 wherein said computer 
program further comprises: 

means for assigning a standard input blood flow and 
standard vessel volume for each vessel compartment; 
5 means for assigning a standard input blood flow, a 

standard capillary volume, and a standard extravascular 
voliune for each organ compartment; 

means for adjusting the standard input blood flow of 
each vessel compartment and each organ compartment as 
10 being proportional to a ratio of cardiac output of the 
patient to a standard cardiac output; 

means for adjusting the standard capillary volume 
and the standard extravascular volume of each organ 
compartment as being proportional to the ratio of blood 
15 volume of the patient to a standard blood volume. 

19. The computer of Claim 18 wherein said computer 
program further comprises means for further subdividing 
each organ into a capillary compartment and an 
extravascular compartment and the step of assigning at 

5 least one differential equation to each organ compartment 
includes assigning at least one differential equation to 
each capillary compartment and each extravascular 
compartment • 

20. The computer of Claim 19 wherein said computer 
progremi further comprises means for converting the 
computed contrast concentration into computed tomography 
enhancement units and means for displaying the 

5 enhancement level of the organ in computed tomography 
enhancements units as a function of time. 
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21. A compu^ted tomography machine comprising: 
a compu-ted -tomography machine; 

a computer having a memory connected to said 
computed tomography machine £or controlling Its 
5 operation ; and 

a computer program In said computer memory, said 
computer program having means for predicting a structure 
specific CT enhancement level In a patient for a given 
patient habitus and a specific contrast Injection 
10 protocol. 

22. The computed tomography machine of Claim 21, 
wherein said computer program comprises: 

means for predicting prior to Implementing an 
Injection protocol an organ enhancement level as a 
5 function of time elapsed after Implementing the Injection 
protocol based on t:he Injection protocol and patlen-t 
specific parameters; 

means for determining an optimum scan start time and 
scan duration based on the predicted enhancement level; 
10 and 

means for performing the scan In accordance with the 
pre-determlned scan start time and scan duration. 

23* A method of determining a set of parameters for 
an Injection protocol for scanning a tissue In a patient: 
using computed tomography, t:he patient having a plurality 
of patient specific parameters, wherein the tissue to be 
5 scanned Is enhanced with an Intravascular ly Injected 
contrast agent, comprising the steps of: 

generating a tissue enhancement function comprising 
a tissue enhancement level for the tissue to be scanned 
as a function of an elapsed time after Injection based on 
10 the patient specific parameters and a specified Injection 
protocol ; and 

determining the set of Injection protocol parameters 
for an optimum scan based on the predicted tissue 
enhancement function. 



wo 97/12550 



PCT/US96/15680 



60 

24. The method of Claim 23 wherein the step of 
determining the set of parameters for an optimum scan 
includes using the predicted tissue enhancement function 
to thereby determine an injection flow rate and a 

5 contrast volume which is predicted to cause the tissue 
enhancement level to exceed a pre- selected threshold 
value • 

25. The method of Claim 24 further comprising the 
step of adjusting the set of parameters until the tissue 
enhancement level is predicted to exceed the threshold 
value* 

26. The method of Claim 25 wherein the set of 
parameters is sequentially adjusted until the predicted 
tissue enhancement function is maintained above the 
threshold value for a time period at least approximately 

5 equal to a specified scan duration • 

27. The method of Claim 26 wherein the step of 
adjusting the set of parameters includes selecting the 
injection flow rate and selecting the injection volume* 

28. The method of Claim 23 wherein the step of 
determining a set of parameters for an optimum scan 
includes determining an optimum scan interval equal to a 
specified scan duration during which l:he predicted tissue 

5 enhancement function is the greatest* 

29. The method of Claim 28 wherein the step of 
determining the optimum scan interval includes 
integrating the predicted enhancement function for 
successive intervals equal to the specified scan duration 

5 and selecting an interval having the greatest integration 
value as the optimum scan interval. 

30. The method of Claim 23, further comprising the 
steps of: 

predicting prior to injecting the contrast agent a 
regional enhancement level as a function of an elapsed 
5 time after injection for a region of interest in the 
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patient based on the patient specific parameters and 
specified injection protocol; 

initiating the injection protocol; 

sequentially measuring an actual enhancement level 
10 in the region of interest at predetermined elapsed times 
after initiating the injection to generate an actual 
regional enhancement function; 

calculating a correction factor based on the 
relationship between the predicted regional enhancement 
15 function and the actual regional enhancement function; 
and 

using the correction factor to calibrate the 
predicted tissue enhancement function. 

31 • The method of Claim 30 wherein at least one of 
the plurality of patient specific parameters has an 
unknown value and wherein the step of predicting a 
regional enhancement function includes the steps of: 
5 providing a mathematical model of a cardiovascular 

system, the model mathematically describing transport of 
the contrast agent through the cardiovascular system; 

inputting a plurality of alternatives for the 
unknown patient parameter into the model to generate a 
10 set of alternative regional enhancement functions with 
the set comprising a member for each alternative; and 

selecting a predicted regional enhancement function 
from the set of alternative regional enhancement 
functions • 

32, The method of Claim 31 wherein the step of 
selecting a predicted regional enhancement function from 
the set of alternative regional enhancement functions 
comprises the steps of: 
5 comparing the actual regional enhancement function 

to the members of the set of alternative regional 
enhancement functions; and 
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selecting from the set of alternative regional 
enhancement functions the member which most closely 
10 resembles the actual enhancement function. 

33, The method of Claim 32 wherein the step of 
selecting the member from the set of alternative regional 
enhancement functions includes the steps of: 

graphing the members of the set of alternative 
5 regional enhancement functions; 

graphing the actual regional enhancement functions- 
integrating the actual regional enhancement graph 

from start to a predetermined time after injection; 

determining which member of the alternative regional 
10 enhancement graphs at the predetermined time has an 

integrated area closest to the integrated area of the 

actual regional enhancement graph at the predetermined 

time. 

34. The method of Claim 31 wherein the tissue to be 
scanned is distinct from the region of interest and 
wherein the step of using the correction factor to 
calibrate the predicted tissue enhancement function 

5 includes the steps of: 

Inputting into the mathematical model the 
alternative chosen for the unknown patient parameter 
which was used to generate the predicted regional 
enhancement function; and 
10 generating a revised predicted tissue enhancement 

function. 

35. The method of Claim 34 wherein the step of 
determining a set of parameters for an optimum scan 
Includes the step of: 

using the revised predicted tissue enhancement 
function to determine an optimum scan delay to cause the 
5 onset of the scanning of the tissue to coincide with a 
threshold enhancement level in the tissue. 

36, The method of Claim 30 wherein the step of 
sequentially measuring an actual enhancement level is 
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performed using low-dose pre- scanning of the region of 
interest. 

37. The method of Claim 31 wherein the region of 
interest is the tissue to be scanned. 

38. The method of Claim 31 wherein said 
mathematical model is embodied in a computer program. 

39. A contrast injector system comprising: 
a contrast injector; 

a computer having a memory connected to said 
contrast injector for controlling its operation; and 
5 a computer program in said computer memory, said 

computer program having means for predicting a structure 
specific CT enhancement level in a patient having a 
specific patient habitus based on a specified injection 
protocol • 

40. The contrast injector of Claim 39, wherein said 
computer program comprises: 

means for predicting prior to implementing the 
injection protocol a tissue enhancement level as a 
5 function of elapsed time based on the injection protocol 
and the specific patient habitus; and 

means for determining an optimum injection flow rate 
and an optimum contrast volume based on the predicted 
tissue enhancement level. 

41. The contrast injector of Claim 40 wherein the 
computer program further comprises: 

means for accepting an input including a revised 
injection rate and a revised contrast volume; 
5 means for revising the predicted tissue enhancement 

level based on the revised injection flow rate and the 
revised contrast vol\ime. 

42. The contrast injector of Claim 41 wherein the 
computer program includes means for Implementing the 
injection protocol using the revised injection flow rate 
and revised contrast volume. 
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43. The contrast Injector of Claim 42 wherein the 
contrast Injector Is In communication with a computed 
tomography machine and wherein the computer program 
Includes means for signaling the computed tomography 

5 machine to perform a scan. 

44. A computed tomography system for scanning a 
tissue In a patient, the patient having a plurality of 
patient specific parameters, wherein the tissue to be 
scanned Is enhanced with an intravascular ly Injected 

5 contrast agent, the system comprising: 
a computed tomography machine; 
an Injector system in communication with the 
computed tomography machine; 

a computer having a memory, the computer in 
10 communication with both the injector system and the 
computed tomography machine for coordinating their 
operation; and 

a computer program in the computer memory, the 
computer program having means for predicting prior to 
15 Injecting the contrast agent a tissue enhancement level 
for the tissue based on the patient specific parameters 
and a selected injection protocol. 

45. The system of Claim 44 wherein the computer 
program includes: 

means for determining an optimum injection protocol 
for injection of the contrast agent; and 
5 means for implementing the optimum injection 

protocol • 

46. The system of Claim 44 further comprising means 
for monitoring actual regional enhancement levels In a 
region of Interest In the patient after initiation of the 
Injection protocol. 

47. The system of Claim 46 wherein the computer 
program includes means for receiving and analyzing the 
regional enhancement levels to thereby confirm the 
predicted tissue enhancement level. 
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48. The sys-tem of Claim 47 wherein the computer 
program Includes means for revising the predicted tissue 
enhancement level based on the regional enhancement 
level . 

49. A computer being prograiraned for predlct:lng a 
structure specific CT enhancement: level In a patient for 
a given patient habitus and a specific contrast Injection 
protocol comprising : 

5 a computer having a memory; 

a computer program In said memory^ the program 
having means for accepting patient specific Information 
and contrast specific Information which affect 
physiological parameters of contrast enhancement: and 
10 means for computing and output ting structure specific 

enhancement levels as a funct:ion of an elapsed t:lme after 
injection. 

50. The computer of Claim 49 wherein t:he comput:er 
program further comprises: 

means for optimizing the injection protocol to 
thereby cause at least a threshold contrast enhancement 
5 level over a predetermined scan duration; and 

means for communicating the optimum injection 
protocol to an Injector. 

51. The computer of Claim 49 wherein the computer 
program Includes; 

means for determining a set of scan paramet:ers for 
an optimum scan based on the computed enhancement levels; 
5 means for communicating the set of scan parameters 

to a CT scan machine. 

52. The computer of Claim 49 wherein the computer 
program Includes means for monitoring actual enhancement 
levels In a region of Interest to thereby confirm the 
computed enhancement levels. 

53. A method for determining an Injection prot:ocol 
for Injecting a contrast agent Into a patient for 
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enhancing a tissue desired t:o be scanned with computed 
tomography, the method comprising the steps of: 
5 generating a plurality of predicted tissue 

enhancement functions based on different patient 
parameters; and 

selecting one of said predicted tissue enhancement 
functions as being close to expected based upon feedback 
10 from a low dose scan of a region of interest. 

54. The method of Claim 53 wherein the selecting 
step Includes the steps of generating from said low dose 
scan an actual tissue enhancement function of said region 
of Interest, Integrating said actual tissue enhancement 

5 function over a predetermined time period, comparing said 
Integrated value with corresponding integrated values for 
each of the predicted tissue enhancement functions and 
selecting as desired the predicted tissue enhancement 
function with the closest Integrated value. 

55. The method of Claim 54 further comprising the 
steps of determining a patient parameter based upon the 
selected predicted tissue enhancement function, and 
determining a predicted tissue enhancement function for 

5 the tissue to be scanned based at least in part on said 
determined patient parameter. 

56. The method of Claim 53 wherein the low dose 
scan is of a tissue different than the tissue desired to 
be scanned. 

57. A method for determining an injection protocol 
for injecting a contrast agent into a patient for 
enhancing a tissue desired to be scanned with computed 
tomography, the method comprising the steps of: 

5 generating a predicted tissue enhancement function; 

selecting a threshold level; and 

adjusting the injection protocol until the tissue 
enhancement function remains greater than the threshold 
level for a time period at least equal to a desired scan 
10 duration. 
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58. The method of Claim 57 further comprising the 
steps of determining an optimal scan start time and end 
time by Integrating the tissue enhancement function for 
successive scan duration Intervals contained within the 

5 portion wherein the threshold is exceeded, and choosing 
as optimal the Integrated interval having the greatest 
value. 

59 . The method of Claim 53 further comprising the 
steps of implementing said method in a computer, and 
performing a CT scan using said determined injection 
protocol. 

60. The method of Claim 57 further comprising the 
steps of implementing said method in a computer, and 
performing a CT scan using said determined injection 
protocol • 
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which vUl afPect physiological 
paraaeters oP contrast enhanceaent 



accept peraeabilitY Factor 
Cany value Proa 2 to 10 



user irpjt 



accept body size option 
(standard or user input) 



standard 



Accept patient inFornatlon 

(gendep.wel^t. height. age.cardloc output index) 

and aodlPY blood voluae and cardiac output 



standard blood volute 
sSQOO wL 

standard cardiac output 
-6S00 aL/atn 



blood voluae (aL) Per aAilt aale 

»33. IMU heightiiO. 725 »( velghtliO. 125 ) 

blood voluae ( aL ) Per adult Peaale 

-33. 161l( helghtllO. 725 )I( veightIlO.425 ) 

cardiac output aL/ain ) 

836.36« helghtllO. 725 )l( weightnO. 425 ) 



I 



ir age >30. cardiac output 
-cardiac output i 
(l-0.0Q8i(age-30)} 



■ultiply cardiac output index: 
(norB8l:l.lov:<l.high:>l ). It 
aay be deterained froi sequential 
CT scanning or the aorta 



I Accept contrast concentration 



Accept injection aethod 
(aano-v8 bi phasic) 



bi-phaslc 



Accept injection (tiael and tiae2) 
and injection (ratal and rate2) 



Bono-phaslc 



Accept total injection tiae 
and injection rate 



return 



FIG . 7 
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Subroutine coip8rtwnt.Mdel : 



starting froa the right heart. 
Pollov the blood riou in the 
body BOdel diagraa 







assign standard input blood flow ( Ol ) 
and vessel voluae ( vi ) ror each eleaent 
structire for e vessel coapartaent ( 1 ) 






assign standard inpu 
capillary voIum CVJ 
voluae (Vj^l) ror eai 
for an organ coapart 


t blood flow (Qj ). 
1. and extravascular 
ch eleaent structure 
aent Cj ) 



VI 



01 



VJ 

VJ*1 



OJ 



adjust the blood Plow and voluae of eech 
coapartaent proportional to the ratio or 
the cardiac output and the blood voluae 
or e patient to those or the standard 



Assign a dirrerentlal equation describing 
contrast trar^port in a vessel coapartaent 
1 1 ):dCl/dtB((Oi-liCl-l )-(OllCl l)/Vi 



Ol-l.Cl-1 



VI 



01. Ci 



Assign e dirrerentlal equation describing 
contrast transport in a caplUarv coapartaent 
(j ):dCj/dt-((Oj-llCj-l HOJICJ »)/vj- 

peraeabl I ity factoriOjw Cj -Cj ♦! )/vj 
and extravasculer coapertaent (j«l):dCj*l/dt 
speraeablllty ractoriOjl(Cj-Cj«l)/yj«i 



Oj-l.CJ-1 



VJ 

vjn 



Qj.cj 



return 



FIG . 8 
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I M I I I I I I I I I 1 I I I I I I I I I I I I I i I I 
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Concentration ( 1 ig/cc ) 



350-1 



300- 



250- 



^200-1 
c 

ilSO 
(J 

c 
o 

c 100- 



50- 



FIG . 9 



- o Aorta 
P Liver 
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FIG. 1 □( b ) 
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FIG . lie a ) 




FIG . lie b ) 
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HO ^250- 
I 200^ 




'I I ' I I I I I I I I I I I I I I I I 
0 30 60 90 120 150180 210240 270 300 
Time ( sec ) 



FIG . 12( a ) 



-,350n 



Aorta 
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' I I I I I I I I I I I I I I I I I I 
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Time ( sec ) 
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500 

= 400- 

5350- 

■K300- 

|250- 

o200- 
t- 



c 

CD 
£ 

C 
Ul 



^ ISO- 
is 100- 

50- 



110 lb8 
ISO lbs 
200 lbs 
250 lbs 



1 I — I — I — I — I — I — I — I — I — (— 1 — I — I — I I I I I I 
0 0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 
Time ( min ) 
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0 0.5 1 1.5 2 2.5 3 3.5 ^ -1.5 5 
Time ( min ) 

FIG. 1 4( a ) 



0.1 




0 0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 
Time ( min ) 



FIG. 1 4C b ) 
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Enhancement (H.u. > 



TIME (MIN. ) AORTIC HEPATIC 



0 


0 


0 


0.13 


15.23 


0.02 


0.25 


187.3 


3. 12 


0.38 


285.65 


15.53 


0 .51 


324 . 22 


34.8 


0.64 


176.9 


52.39 


0.77 


135.62 


60.29 


0.91 


127.09 


62.9 


1.04 


127.15 


63.57 


1.17 


122.31 


63.57 


1.3 


114.23 


63 


1.43 


108.87 


61.84 


1.56 


106.65 


60.38 


1.69 


104.54 


58.93 


1.82 


100.87 


57.47 


1.95 


96.48 


56 


2.07 


92.45 


54.45 


2. 2 


89.1 


52.81 


2.34 


86.21 


51.15 


2.47 


83.57 


49.55 


2.59 


81.09 


48.04 


2. 73 


78.65 


46.51 


2.85 


76.5 


45. 11 


2.98 


74.49 


43.79 


3. 11 


72.55 


42.51 


3. 23 


70. 72 


41.29 


^ . J / 


bo . 93 


40. 1 


3.49 


67.25 


38.97 


3.63 


65.53 


37.87 


3. 75 


64.17 


36.87 


3.88 


62.78 


35.93 


4.01 


61.45 


35.01 


4. 14 


60.21 


34. 16 


4.57 


58.99 


33. 3^ 


4.4 


57.83 


32.52 


4.52 


56.75 


31.78 


4.65 


55.75 


31. 1 


4 .78 


54.79 


30. 44 


4.9 


53.85 


29.79 


5 


S3. 18 


29.34 
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FIG . 17 



wo 97/12550 



PCT/US94in5680 



ENHANCEMENT 



14/20 

FIG . 18 

•152 




THRESHOLD 



TIME 



ENHANCEMENT 



FIG. 19 




•MAX. AREA UNDER CURVE ABOVE 
THRESHOLD GIVEN INTERVAL ( D-C ) 

THRESHOLD 



TIME 



ENHANCEMENT 
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PATIENT 



A6E 

GENOER 
HEIGHT 
WEIGHT 

CAROIAC OUTPUT 



(GOOD. FAIR. POOR) 



INJECTION 



UNI- vs. BIPHASIC 
VOLUME 

CONCENTRATION 
INJECTION RATE 



UNI 
120al 

300 igl/il 
3al/8ec 



CT SCAN 



EXAN TYPE (LIVER. 
SCAN DURATION 



VASCULAR) 



LIVER 
30 seconds 



OUTPUT PARAMETERS; 

PREOICTEO CONTRAST ENHANCEMENT CURVE 

EIMANCEMENT ADEQUACY (y/n) 

SCAN DURATION ABOVE ADEQUACY THRESHOLD 
IF NO. INCREASE VOLUME. OR RATE 
IF YES. DECREASE VOLUME. OR RATE 



OPTIMAL SCAN OEUY 

CARDIAC OUTPUT - lOOX 



TIME 


PREOICTEO 


SUBTRACT 


COMPUTE AREA 


OPTIMAL 




HEPATIC 


THRESHXO OF 


UNDER CURVE (AUC) 


TEMPORAL 




ENHANCEMENT 


ADEQUACY 


OVER SCAN DURATION 




LEVEL 


(SO HU) 


( 30 SECONDS ) 


UINOOU 


SECONDS 


HU 


HU 


HUoSEC 




0 


0 


-50 






5 


0 


-30 






10 


0 


-50 






IS 


1 


-49 






20 




-46 






23 


10 


•40 






30 


17 


-33 






35 


25 


•25 






iO 


34 


-16 


9.6 




15 


43 


-7 


129.0 




SO 


49 


-1 


205.6 




55 


53 


3 


246.6 




EO 


56 


6 


263.8 




B5 


56 


8 


263.7 


START SCAN 


70 


SB 


B 


252.6 




15 


SB 


B 


233.9 




BO 


SB 


8 


209.8 




BS 


SB 


8 


181.9 




BO 


S7 


7 


151.4 


END SCAN 


95 


56 


6 


119.4 




100 


SS 


5 


86.0 




105 


54 


4 


SI. 5 





SCAN IS STARTED AT MAXIMUM AUC (263.8 HU*SEC) OVER THE SCAN OURATIOM 
(BOLD " SCAN DURATION) 

FIG . 20 
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START 



INPUT: CONTRAST ENHANCE CURVES 

DESIRED ENHANCEICNT THRESHOLI^ 
SCAN OtAATION 



RE-RUN FOR NEU CONTRAST 
ENHANCENENT CURVES 



LOWER THE THRESHOLD 




I 



IF C VOLUME m\ 
RATE < HUU 



IF (DURATION ABOVE THRESHOLD 
> = SCAN DURATION) 



t VOLUME OR 

t INJECTION RATE 



I 



IF (DURATION ABOVE THRESHOLD 
" SCAN DURATION! 



IF (VOLUHE OR 
RATE > HIN. I 



USER OPTION ( 1 ): 
i VOLUHE OR 
i- INJECTION RATE 



USER OPTION ( 2 h 
EVALUATE OPTIMAL 
SCAM DURATION 



CALCULATE AREAS UNDER THE ENHANCEMENT CURVE 
(AUCl VITH A FIXED INTERVAL OF SCAN DURATION. 
STARTING FROH THE INITIAL TO THE END-POINT OF 
THE DURATION ABOVE THRESHOLD 



DURATION ABOVE THRESHOLD^ 
DPTIIIAL TEiraUL VINOM 



SCAN DURATION PROVIDING THE MAX. AUC: 
OPTIHAL TEHPQRAL VINOOV 



RETURN 



FIG. 21 
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TIME ( min. ) 



FIG . 22A 




O.S 



1.5 2 2.5 3 3.5 4 
TIME (min. ) 



4.5 



FIG . 22B 
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CALCULATION OF THE TEMWAL WINDOW FOR OPTINAL ENHANCEMENT 



CARDIAC OUTPUT - 7SZ 



TINE 



PREDICTED 
HEPATIC 
ENHANCENENT 
LEVEL 



SUBTRACT 
THRESHOLD OF 
ADEQUACY 
tSOHU) 



CONPUTE AREA 
UNDER CURVE (AUC) 
OVER SCAN DURATION 

(30 SECONDS) 



OPTIHAL 
TENPORAL 

WINDOW 



SECONDS 



0 
5 
10 
15 
20 
25 
30 
35 
10 
15 
SO 
55 
60 
65 
70 
75 
BO 
85 
90 
95 
100 
105 
110 
115 
120 
125 
130 
135 
HO 
115 
ISO 
155 
160 
165 
170 



NU 



0 
0 
0 
0 
2 
1 
9 
16 
23 
30 
38 
15 
49 
53 
55 
57 



59 



SB 
SB 

57 

56 

56 

55 

54 

54 

53 

52 

52 

51 

SO 

49 



HU 



-50.0 
-50.0 
-50.0 
-49.7 
-48.5 
-45.6 
-40.7 
-34.3 
-27.4 
-19.8 
-11.6 
-5.2 
-0.5 
2.8 
5.1 
6.6 
7.6 
8.2 
8.5 
8.5 
8.4 
8.0 

7.6 
7.1 
6.4 
5.8 
5.1 
4.4 
3.7 
3.0 
2.3 
1.6 
0.9 
0.1 
-0.7 



HUoSEC 



23.8 
122.8 
191.0 
236.3 
264.5 
279.2 
284.0 
281.3 
272.6 
259.3 
242.1 
222.3 
200.7 
177,8 
154.2 
130.1 
105.7 
80.6 
54.8 



START SCAN 



END SCAN 



SCAN IS STARTED AT NAXINUN AUC C 284.0 HU»SECI OVER THE SCAN DURATION 
ISOLD - SCAN DURATION) 



FIG . 25 
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SUBROUTINE A0JUSTIN6 ONSET OF SCANNING 
(INOEPEMKNT INJECTOR SYSTEH) 



INPUT: PKCIOTED UITH STANDARD CAflOIAC OUTPUT 
AUCl (AAEA IMER EARLY AORTIC CURVE) 
tl (SCAN DURATION TO ACHIEVE AUCl ) 
tIKPREOICTEO ONSET OF HEPATIC SCAN) 



FOR DIFFERENT LEVELS OF CARDIAC OUTPUT. 
CALCULATE A SET OF tl'S YIELDING IDENTICAL 
AUCl AND CORRESPONDING tll'S 



CONSTRUCT AN INTERPOLATING LOOK-UP TABLE 
FOR tll'S AT VARIOUS tl'S 



HEASUREyENT FIUM UV-OOSE PflESCAN CT 



CALCULATE CUHULATIVE AUCl (Ct-tl ) 
AS SCAN PHQCEEOS ( INCREMENTAL At I 



CiF (AUCl 8Aucr) -iL 

Y 



I 



t2l (UPDATED ONSET OF HEPATIC SCAN) 
IS DETERMINED FROM THE LOOK-UP TABLE 



RETURN 



FIG . 26 
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